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SECTION  1 


INTRODUCTION  AND  SUMMARY 


1.1  INTRODUCTION 

A  requirement  of  underground  nuclear  testing  is  that  radioactive 
gases  be  prevented  from  entering  the  atmosphere.  In  general,  this 
requirement  will  be  met  if  experimental  tunnels  leading  from  the  nuclear 
device  cavity  are  stemmed  successfully  and  if  the  cavity  gases  are 
contained  by  the  adjacent  surrounding  medium.  The  residual  stress  field 
generated  in  this  surrounding  medium  by  the  explosion  probably  aids 
containment.  Although  containment  has  been  achieved  for  many  years  and 
stemming  has  generally  been  successful  in  recent  nuclear  tests,  reliability 
is  uncertain  and  planned  tests  still  require  extensive  containment  evaluation 

Over  the  years,  the  DNA  stemming  and  containment  (SAC)  program  has 

consisted  of  five  main  parts:  code  development  for  ground  motion,  tunnel 

closure,  and  grout  flow  calculations;  material  properties  determination; 

laboratory  investigations;  scaled  higli  explosive  tests;  and  field  diagnostics 

Over  the  past  several  years  SRI  has  been  conducting  the  laboratory 
1-6 

investigations . 

Laboratory  investigations  during  the  last  year  have  focused  mainly 
on  containment.  One  purpose  of  the  experimental  program  lias  been  to 

compare  the  containment  capability  of  cavities  generated  under  a  variety 

.  .  7-9 

of  conditions.  Another  purpose  is  to  validate  calculations  of  the 

residual  stress  field  as  a  major  containment  feature  of  underground  nuclear 

tests.  The  data  generated  in  the  laboratory  are  suitable  for  correlation 

with  the  predictions  of  existing  codes.  The  experimental  technique  has 

also  been  used  to  determine  the  influence  of  major  geological  and  test 

site  features.  Although  our  basic  containment  experiment  may  not  be 

regarded  as  a  small-scale  version  of  an  explosively  simulated  underground 

nuclear  test,  simulation  is  close  enough  to  provide  the  correct  important 

mechanisms  required  in  a  study  of  containment. 
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In  our  basic  containment  experiment,  a  small  spherical  charge  is 
cast  at  the  center  of  a  sphere  of  rock-matching  grout  (RMG)  with  properties 
similar  to  those  of  Nevada  Test  Site  tuff.  The  RMG  sphere  is  hydraulically 
pressurized  to  represent  overburden.  The  explosive  is  detonated,  and  dyed 
fluid  is  immediately  pumped  into  the  cavity  at  a  constant  flow  rate  until 
fracture  occurs  and  a  steady  radial  pressure  gradient  is  established. 

(In  an  alternative  version  of  this  test,  residual  cavity  gases  are  vented 
before  the  start  of  hydrofracture.)  The  internal  pressure  is  removed, 
and  then  the  external  pressure  is  removed.  The  cracked  sphere  is  drained 
and  tapped  into  two  parts  with  a  chisel,  and  the  exposed  fracture  surfaces 
are  photographed.  The  effect  of  the  explosively  generated  residual  stress 
field  on  containment  is  assessed  by  conducting  a  separate  hydrofracture 
test  on  a  sphere  with  an  unexploded  cavity  equal  in  size  to  the  exploded 
cavity. 

Exploded  and  unexploded  cavity  hydrofracture  tests  were  performed  on 

RMG  spheres  of  12-inch  (30.48-cm)  diameter  to  establish  reproducibility 

of  the  basic  experimental  technique.  In  exploded  cavity  tests,  a  3/8-gram 

charge  of  PETN  generated  a  3/4-inch  (1.90~cm)  diameter  cavity,  which  was 

unvented.  In  unexploded  cavity  tests,  a  smooth  and  unlined  cavity  having 

the  same  diameter  as  the  explosively  formed  cavity  was  cast.  The  fixed 

parameters  were  overburden  pressure  [1000  psi  (6.895  MPa)],  rate  of  fluid 

3 

flow  into  the  cavity  (4.26  cm  /min),  and  viscosity  of  the  hydrofracture 
fluid  (1  centipoi.se).  Surface  gages  monitored  fracture  and  fluid  arrival. 
Embedded  gages  monitored  fluid  arrival.  Results  of  the  reproducibility 
tests  provide  comparison  data  for  assessing  the  effects  on  containment  of 
the  following  features: 

•  Overburden.  External  pressures  of  0,  1000, 
and  2000  psi  (0,  6.895,  and  13.79  MPa)  were 
applied  to  unexploded  cavity  spheres. 

•  Viscosity.  Unexploded  and  exploded  cavity  spheres  were 
hydrof ractured  with  water  (1  centipoise  viscosity)  and 
glycerol  (660  centipoise  viscosity). 

£ 

Pentaerythr itol  tetranitratc  (C  H  0  ,N  ) . 
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•  Material  property.  Rock-match ing  grout  (RMG  2C4), 
low  density  rock-matching  grout  (LD  2C4),  and 
granite  simulant  ((.’S3)  spheres  were  hydrof rartured . 

•  Cavity  lining.  Hydrof raeture  tests  were  performed  on 
unexploded  cavities  lined  with  an  impermeable  membrane. 

•  Charge  holder  thickness.  l.ucite  charge  holders  with 
lb-mil  (0.406-mm)  and  26-mil  (0.660-mm)  wall  thicknesses 
wire  cast  in  exploded  cavity  spheres. 

•  I’neoupled  charge.  A  3/16-inch  (0.476-cm)  air  gap 
between  charge  and  surrounding  grout  was  provided 
in  exploded  cavity  tests. 

•  Geometry.  exploded  cavity  cylinders  12  inches  (30.48  cm) 
in  diameter  and  12  inches  (30.48  cm)  long  were 

hydrof rac tured . 

•  Material  interface.  Planes  of  discontinuity  in  material 
properties  were  formed  in  exploded  cavity  spheres  by 
bonding  rock-matching  grout  (RMG  2C4)  to  low  density 
rock-matching  grout  (LD  2C4)  and  by  bonding  RMG  2C4 

to  granite  simulant  (GS3). 

•  Geological  and  test  site  features.  Radial  tunnels, 
nearby  tunnels,  and  fault  planes  were  east  in  exploded 
cavity  spheres. 

•  Stress  relaxation.  Variations  were  permitted  in  the 
decay  of  residual  cavity  gas  pressure  in  unvented 
exploded  cavities. 

In  addition,  the  reproducibility  of  the  3/8-gram  explosive  source 
was  assessed  by  monitoring  the  pressure  pulse  generated  in  the  overburden 
fluid  during  exploded  cavity  tests.  A  fracture  mechanics  analysis  was 
applied  to  unexploded  cavity  tests  to  explain  the  dependence  of  fracture 
initiation  pressure  on  hydrof raeture  fluid  viscosity  and  overburden.  A 
porous  flow  analysis  was  applied  to  exploded  cavity  tests  to  assess  the 
rate  of  decay  of  the  explosively  generated  residual  stress  field.  Particle 
velocity  measurements  obtained  in  an  independent  laboratory  program'11 
provided  additional  data  for  verifying  the  calculations. 

Tin1  tests  performed  during  the  current  phase  . . r  the  laboratory  program 
are  tabulated  according  to  category  in  Section  3.1.  In  general,  a  minimum 
of  three  tests  in  each  category  is  desirable  to  demonstrate  reproducibility 
or  to  show  a  trend  developing.  However,  the  extensive  variety  of  experi¬ 
mental  configurations  limited  the  number  of  tests  in  a  few  categories  to 


one  or  two.  Section  2  describes  the  most  recent  experimental  developments, 
and  Section  3  presents  and  discusses  the  hydrofracture  results.  Appendix  A 
describes  the  application  of  a  fracture  mechanics  analysis  to  hydrofracture. 
Appendix  B  examines  relaxation  of  a  residual  stress  field  in  terms  of 
porous  flow  theory.  Appendix  C  describes  material  properties  for  rock¬ 
matching  grout  RMG  2C4,  low  density  rock-matching  grout  LD  2C4,  and  granite 
simulant  GS3.  Appendix  D  describes  the  sequence  of  events  in  unvented 
exploded  cavity  tests. 


1 . 2  SUMMARY 

The  principal  findings  of  the  above  investigations  were  as  follows: 

•  Reproducibility  of  the  unexploded  cavity  hydrofracture 
records  is  very  good  (Figure  3.1). 

•  Reproducibility  of  the  explosive  source  is  good  for  both, 
charge  holder  wall  thicknesses  tested.  The  wall  thickness 
variation  has  negligible  effect  on  hydrofracture. 

•  Reproducibility  of  unvented  exploded  cavities,  measured  in 
terms  of  fracture  initiation  pressures,  is  good  for  water 
and  glycerol  hydrofractures.  The  pressure  records  depend 
on  the  residual  cavity  gas  pressure  at  the  start  of 
hydrofracture. 

•  Fracture  initiation  pressures  for  water  hydrofracture 
of  unvented  exploded  cavities  are  nearly  twice  the 
corresponding  values  found  in  unexploded  cavity  tests. 

•  Hydrofracture  pressures  for  vented  exploded  cavity 
spheres  span  the  range  between  those  for  unexploded 
cavity  tests  and  unvented  exploded  cavity  tests. 

•  The  slope  of  the  fracture  initiation  pressure  versus 
overburden  curve  is  unity  for  water  and  glycerol 
hydrofractures  of  unexploded  cavity  rock-matciiing 
grout  (RMG  2G4)  spheres. 

•  Increased  hydrofracture  fluid  viscosity  results  in 
higher  fracture  pressures  for  unexploded  cavity  rock¬ 
matching  grout  (RMG  2C4)  and  granite  simulant 

(GS3)  spheres. 

•  Fining  an  unexploded  cavity  with  an  impermeable  membrane 
results  in  a  higher  fracture  initiation  pressure  and 
initially  stable  crack  growth. 

•  Higher  fracture  initiation  pressures  are  developed  in 
unexploded  cavities  surrounded  by  stronger  materials 
such  as  GS3 . 
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•  Air  voids  in  low  density  grout  produee  l.irgor  exploded 
cavities,  less  residual  stress,  and  lower  hydro! rucLure 
pressures . 

•  Uncoupling  the  explosive  source  from  the  surrounding 
medium  has  negligible  effect  on  hvdrof racture  pressures. 

•  Cylindrical  exploded  cavity  specimens  yield  pressure 
records  comparable  to  those  obtained  with  spheres. 

•  A  plane  interface  connecting  RfIC  2C4  and  1.1)  2C4  has 
negligible  effect  on  hvdrofracture.  The  effect  of  a 
plane  interface  connecting  R>1C  2C4  and  ('.S3  is  determined 
by  the  strength  of  the  bond  between  the  materials. 

•  Hollow  radial  tunnels  reduce  hydrofracture  pressure 
because  cavity  gases  vent  to  the  tunnel  which  provides 
an  axis  for  a  fracture  to  initiate. 

•  Explosive  products  have  been  contained  in  tests  with  a 
tunnel  passing  within  one  cavity  radius  of  an 
exploded  cavity. 

•  Containment  has  been  achieved  in  tests  with  a  fault 
plane  passing  within  three  cavity  radii  of  an 
exploded  cavity. 

in  summary,  an  explosively  formed  residual  stress  field  enhances  the 
containment  capability  of  a  cavity  by  increasing  hvdrofracture  pressures 
in  unvented  tests  and  by  stabilizing  initial  crack  growth. 

Investigations  during  the  next  phase  of  the  containment  program  will 
concentrate  on  the  influence  of  pore  water  migration  and  creep  on  the 
residual  stress  field.  The  program  will  include: 

•  Characterizing  the  residual  stress  field  surrounding 
an  exploded  cavity  by  means  of  stress,  strain,  and 
particle  velocity  measurements. 

•  Decreasing  the  time  from  charge  detonation  to  fracture 
initiation  in  (invented  exploded  cavity  tests. 

•  Determining  the  effects  on  containment  of  such  test  site 
features  as  nearby  chimneys  and  asymmetric  overburden  loads. 


SECTION  2 


EXPERIMENTAL  DEVELOPMENT 


2 . 1  CONCEPT 

The  experiment  shown  schematically  in  Figure  2.1  was  devised'*  to 
simulate'  in  the  laboratory  the  conditions  in  an  underground  nuclear  test. 

A  small  spherical  explosive  charge  representing  the  device  is  sealed  in 
a  thin  I.ucite  container  and  cast  in  a  much  larger  sphere  of  rock-matching 
group  (R11G)  with  properties  similar  to  those  of  Nevada  Test  Site  tuff 
[Figure  2.1(a)).  The  grout  sphere  is  cured  and  then  hydraulically  pres¬ 
surized  to  represent  overburden  [Figure  2.1(b)].  The  explosive  is  detonated 
and  the  residual  gas  pressure  is  measured  [Figure  2. 1(c) 1.  Fluid  is 
immediately  pumped  into  the  cavity  at  a  constant  flow  rate  until  fracture 
occurs  and  a  steady  flow  develops  along  the  fracture  surface  [Figure  2.1 
(d) ] .  An  alternative  version  of  this  experiment  was  developed  to  allow 
hydrofracture  following  the  venting  of  cavity  gases. 

The  hydrofracture  portion  of  the  experiment  is  included  to  determine 
the  effect  of  the  residual  stress  field  by  comparing  the  cavity  pressures 
required  to  crack  the  RMG  sphere  with  and  without  residual  stresses.  In 
the  experiments  without  a  residual  stress  field,  spherical  cavities  are 
cast  in  the  grout  sphere;  these  unexploded  cavities  are  the  same  size  as 
the  exploded  cavities. 

2.2  EXPERIMENTAL  TECHNIQUES 

Figure  2.2  shows  the  experimental  apparatus.  A  1 2-inch--d iameter  RMG 
sphere  is  shown  inside  a  steel  vessel  containing  water  that  can  be  pres¬ 
surized  to  the  desired  overburden.  The  sphere  is  suspended  from  the  lid 
by  a  steel  tube  cast  in  the  grout.  In  a  vented  exploded  cavity  test, 
this  tube  provides  a  means  for  positioning  the  charge,  potting  in  the 
detonator  cables,  and  drilling  into  the  cavity  after  detonation.  The 
water  in  the  vessel  is  maintained  at  a  constant  pressure  throughout 


V* 


(b)  OVERBURDEN 


MA-3702-103B 

FIGURE  2.1  SEQUENCE  OF  OPERATIONS  IN  CONTAINMENT  EXPERIMENTS 
(Unvented  Exploded  Cavity} 
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the  test  by  incorporating  a  high-pressure  gas  reservoir  and  valve  in 
the  water  supply  line. 


The  constant  flow  rate  system  shown  schematically  in  Figure  2.3 
conforms  with  standard  hydrofracture  practice.  The  specifications  of 
this  system  are  as  follows: 

Motor:  Dayton  Gear  motor  Model  5K933, 

1/4  hp,  6  rpm,  600  in-lb  torque 

Pump:  High  Pressure  Equipment  Co.  Model  87-6-5 

60-cm3  capacity,  5000-psi  maximum  pressure 
4.26-cm^/min  flow  rate  (0 . 71-cm^/revolution) 

System:  Fluid:  dyed  conductive  water  or  dyed  glycerin 

Volume  (excluding  pump):  31.33  cm^ 

The  motor  shaft  rotates  at  a  constant  angular  velocity  of  6  rpm.  A 

slotted  tubular  coupling  allows  an  extension  of  the  shaft  to  translate 

as  well  as  rotate.  The  shaft  extension  is  threaded  and  rotates  in  a 

fixed  threaded  bearing,  resulting  in  a  constant  velocity  translation 

of  the  pump  piston.  Hence,  fluid  is  driven  at  a  constant  rate  of  4.26 

cmVmin  into  the  cavity  supply  line  (Figure  2.2).  Fluid  flow  is  recorded 

automatically  with  pen  and  chart  by  measuring  the  voltage  change  across 

an  angular  potentiometer  geared  to  the  rotating  shaft;  a  constant  chart 

speed  introduces  time,  and  calibration  test  results  relate  flow  to 

voltage  change.  Similarly,  pressure  is  recorded  by  measuring  the  voltage 

change  across  a  calibrated  resistance  pressure  transducer.  Faster  or 

slower  flow  rates  are  possible  through  a  modification  of  the  gear  system. 

Compliance  of  the  hydrofracture  system  is  measured  by  performing 
a  pressure  test  in  which  the  cavity  is  replaced  by  a  rigid  vessel  of 
equivalent  volume.  The  pressure-volume  curve  for  this  test  is  then 
used  to  compensate  for  the  effects  of  system  compliance  on  tiie  hvdro- 
fracture  records. 

Details  of  the  experimental  techniques  for  unvented,  vented,  and 
unexploded  cavity  tests  are  given  below. 


—  -  -■  - - - 
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Unvented  Cavity  Tests 


The  basic  configuration  for  unvented  exploded  cavity  hydrofracture 

tests  is  show  in  Figure  2.4.  Figure  2.5  shows  the  initial  charge  design, 

which  consists  of  3/8-gram  of  PETN  in  a  3/8-inch-OD  (9.52-mm)  Lucite  case 

with  a  wall  thickness  of  16  mils  (0.406  mm).  A  constant  explosive  density 

3 

of  1  S/  cm  is  desired  for  reproducibility;  achievable  machining  accuracy 

results  in  slight  variation  of  PETN  weight  from  charge  to  charge.  Hence 

3/8  gram  is  the  nominal  weight  of  the  explosive.  The  charges  are  assembled 

3 

by  pressing  PETN  to  a  density  of  1  g/cm  into  a  pair  of  mating  Lucite 

hemispheres.  The  bridgewire  assembly  is  then  positioned  as  shown  in 

Figure  2.5,  using  a  notch  in  one  of  the  hemispheres  as  a  guide,  and  the 

* 

hemispheres  are  snapped  together  and  sealed  with  Homalite.  Figure  2.6 
shows  an  improved  charge  design  in  which  the  wall  tiiickness  of  the  Lucite 
charge  holder  was  increased  to  26  mils  (0.660  mm)  to  ensure  stability 
under  anticipated  external  loading.  A  wide-angle  ferrule  is  attached 
near  the  end  of  the  access  tube  to  compress  the  grout  around  the  tube  as 
the  cavity  expands  and  prevent  a  leak  path  to  the  overburden  from  develop¬ 
ing.  Tlie.  charge  is  positioned  by  gluing  the  charge  holder  to  the  steel 
ball  that  serves  as  a  valve  during  the  test;  the  lead  wires  are  brought 
out  of  the  sphere  along  the  plane  of  the  equator. 

Cavity  gas  pressure  is  measured  and  hydro  fracture  performed  by 
following  the  sequence  of  steps  shown  in  Figure  2.7.  Before  charge 
detonation,  the  entire  system  is  filled  with  hydrofracture  fluid.  The 
steel  ball  attached  to  a  valve  stem  seals  the  end  of  the  access  tube 
[Figure  2.7(a)].  Charge  detonation  expands  the  cavity  past  the  end  of 
the  tube  [Figure  2.7(b)],  and  when  the  ball  valve  is  lowered  [Figure  2.7(c)], 
cavity  gases  pressurize  the  hydrofracture  fluid.  Pumping  may  begin 
immediately.  Cavity  gas  and  tiydrof racture  pressures  are  measured  by  a 
pressure  transducer  located  in  the  fluid  supply  line.  Ideally,  charge 
detonation,  the  start  of  pumping,  and  opening  of  the  access  tube  occur 
simultaneously.  In  practice,  however,  initially  high  cavity  gas  pressure 
delays  opening  of  tiie  access  tube  for  approximately  4  seconds,  and  cavity 
pressures  are  not  measured  during  this  time.  File  process  is  described 
schematically  in  Appendix  D. 

•k 
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FIGURE  2.5  EXPLOSIVE  CHARGE 
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FIGURE  2.6  IMPROVED  EXPLOSIVE  CHARGE 
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FIGURE  2.7  SEQUENCE  OF  OPERATIONS  FOR  UNVENTED  EXPLODED  CAVITY  TESTS 
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WATER 


(c)  START  OF  HYDROFRACTURE  PROCESS 

MA-5958-97 


FIGURE  2.7  SEQUENCE  OF  OPERATIONS  FOR  UNVENTED  EXPLODED  CAVITY 
TESTS  (Concluded) 


Accuracy  of  the  cavity  gas  pressure  measurement  was  assessed  by 
performing  a  calibration  test.  In  this  test  the  exploded  cavity  and  gases 
are  replaced  by  a  rigid  vessel  containing  an  equivalent  volume  of  nitrogen 
at  a  known  pressure.  When  the  ball  valve  is  lowered,  the  compliance  of 
the  hydrofracture  system  allows  the  nitrogen  to  expand  and  drop  in  pressure. 
To  eliminate  the  pressure  drop  due  to  compliance,  the  hydrofracture 
system  is  pressurized  to  a  nominal  level  of  300  psi  (2.068  MPa).  This 
prepressurizing  technique  was  incorporated  into  the  (invented  exploded 
cavity  tests. 

Vented  Cavi ty  Tes_ts 

The  basic  configuration  for  vented  exploded  cavity  hvdrofracture 
tests  is  shown  in  Figure  2,8.  The  charge  is  positioned  by  drawing  the 
lead  wire's  through  the  access  tube  and  filling  the  space  between  charge 
holder  and  tube  witli  epoxy.  The  access  tube  is  filled  witii  epoxy.  Venting 
occurs  when  the  tube  is  drilled  out  following  charge  detonation.  The 
drilling  operation  results  in  a  20-minute  delay  between  charge  detonation 
and  hydrofracture.  The  hydrofracture  pressures  are  measured  by  a  pressure 
transducer  located  in  the  fluid  supply  line. 

Uno x_p  1  od  ed  Cavity  T os  t  s 

The  basic  configuration  for  unexpLoded  cavity  hvdrofracture  tests 
is  shown  in  Figure  2.9.  The  method  of  forming  the  central  cavity  is 
determined  by  the  choice  of  a  lined  or  unlined  cavity. 

A  lined  cavity  is  formed  by  first  casting  a  balloon-shaped  rubber 
membrane  with  an  outside  diameter  equal  to  the  desired  cavity  diameter. 

The  neck  of  the  membrane  is  then  stretched  over  the  end  of  the  access 
tube  and  held  in  place  by  means  of  a  Teflon  ferrule  and  epoxy  seal.  The 
membrane  is  filled  witli  water  to  preserve  its  spherical  shape  while  the 
grout  is  poured  and  cured.  After  the  grout  is  cured  and  water  removed 
from  the  cavity,  the  membrane  remains  in  place  to  act  as  a  cavity  lining 
during  hvdrofracture . 


CONSTANT  FLUID 
FLOW  R^TE 


MA-3702-122A 

FIGURE  2.8  BASIC  CONFIGURATION  FOR  VENTED  EXPLODED  CAVITY  TESTS 
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An  unlined  cavity  is  formed  by  first  inserting  a  rubber  membrane 
(balloon)  through  the  access  tube  and  expanding  the  tip  with  water  to 
the  desired  diameter.  A  spherical  shape  is  maintained  by  means  of  a 
wire  cage,  which  is  temporarily  fastened  to  the  end  of  the  tube.  After 
48  hours,  the  membrane  has  assumed  a  permanent  shape  and  the  cage  is 
removed.  The  membrane  and  tube  are  then  positioned  in  a  Lucite  mold  and 
the  grout  is  cast  around  them.  After  the  grout  has  cured,  the  membrane 
is  removed,  leaving  a  smooth  unlined  cavity. 

2.3  INSTRUMENTATION 

The  crack  growth  and  fluid  motion  associated  with  hydrofracture  and 
the  particle  velocity  associated  with  charge  detonation  were  monitored  by 
the  following  sensors: 

•  Surface  fracture  gage  (Figure  2.10) 

•  Surface  fluid  arrival  gage  (Figure  2.10) 

•  Embedded  fluid  arrival  gage. 

In  addition,  visual  observations  of  surface  cracking  and  fluid  arrival 
were  made  for  unexploded  cavity  tests  conducted  without  overburden. 
Finally,  relevant  particle  velocity  records,  which  were  generated  during 
an  independent  laboratory  study, ^  were  incorporated  into  these  results. 

S u r  f  a c  e  Frac  t  ure_G age 

Surface  fracture  of  a  grout  sphere  is  detected  by  means  of  a  1 /lb- 
inch-wide  (1.59-mm)  stripe  of  electrically  conductive  silver-based  paint. 
Two  copper  tabs  are  first  embedded  in  the  surface  of  the  sphere  during 
casting.  After  the  grout  is  cured,  the  paint  stripe  connects  the  tabs  to 
cover  a  desired  region.  During  hydrofracture  the  resistance  from  tab  to 
tab  is  monitored.  Visual  inspection  of  the  surface  of  the  sphere  during 
hydrofracture  has  shown  that  a  hairline  crack  in  the  grout  is  sufficient 
to  break  the  paint  and  produce  a  noticeable  resistance  change  in  the  gage. 
Typical  regions  of  surface  covered  by  the  gage  are  shown  in  Figure  2.10. 
Initial  gage  resistance,  which  is  a  function  of  length,  is  100  to  200  ohms 
for  the  gage  shown.  Final  resistance,  which  is  a  measure  of  electrical 
conductivity  through  the  grout,  is  40,000  ohms. 
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FIGURE  2.10  SURFACE  CRACK  AND  FLUID  ARRIVAL  GAGES 


Surface  Fluid  Arrival  Gage 

Arrival  of  hydrofracture  fluid  at  the  surface  of  a  sphere  is  detected 
by  means  of  a  20-mil  (0.51-mm)  copper  wire  taped  to  the  grout.  The 
resistance  between  access  tube  and  gage  is  monitored  during  a  test. 
Conductive  hydrofracture  fluid  flows  from  the  cavity ,  along  a  crack,  and 
to  the  gage,  resulting  in  a  decrease  in  the  measured  resistance.  The 
portion  of  the  sphere  in  contact  with  the  gage  is  coated  with  a  brittle, 
electrically  insulating  material  to  provide  a  high  initial  resistance  from 
access  tube  to  copper  wire.  With  this  technique,  the  resistance  decreases 
from  1500  to  700  ohms  as  fluid  reaches  the  surface.  Each  copper  wire  may 
be  positioned  adjacent  to  a  surface  paint  stripe,  as  shown  in  Figure  2.10, 
to  allow  for  a  correlation  of  surface  fracture  and  fluid  arrival.  The 
conductive  hydrofracture  fluid  is  a  mixture  of  dyed  water  and  copper 
sulfate. 

Embedded  Fluid  Arrival  Gage 

Arrival  of  hydrofracture  fluid  at  some  interior  point  of  a  grout 
sphere  is  detected  by  means  of  embedded  2-mil  (0.051-mm)  copper  wire. 

The  technique  is  similar  to  that  of  the  surface  fluid  arrival  gage 
described  above.  Conductive  hydrofracture  fluid  flows  along  a  crack 
and  changes  the  resistance  between  the  access  tube  and  the  wire.  Hv 
adjusting  the  configuration  and  orientation  of  the  wire,  we  can  detect 
fluid  arrival  at  any  point  in  the  sphere. 

Embedded  Particle  Velocity  Gage 

Particle  velocity  records  associated  with  charge  detonation  in 
rock-matching  grout  are  shown  in  Figure  2.11.  The  results,  which  were 
taken  from  a  recent  laboratory  study  of  stress  waves  crossing  fault 
planes,  ^  provide  additional  data  for  verifying  the  calculations.  Radial 
velocity  4  cm  from  the  center  of  the  charge  was  measured.  An  average 
profile  based  on  the  six  records  in  Figure  2.11  yields  a  time  of  arrival 
of  15  psec,  a  peak  particle  velocity  of  11.0  m/sec,  and  a  reversal  of 
motion  at  41  psec. 


SO 


Overburden  Pressure  1000  psi  (6.895  MPa) 
Material  RMG  2C4 


FIGURE  2.11  PARTICLE  VELOCITY  ASSOCIATED  WITH  CHARGE  DETONATION  IN  EXPLODED  CAVITY  TESTS 


2.4  CHARGE  REPRODUCIBILITY 


The  reproducibility  of  the  explosive  source  is  a  basic  requirement 
of  exploded  cavity  tests  and  has  been  demonstrated  in  the  past.^  To 
continually  monitor  charge  reproducibility,  we  used  quartz  gages  to 
measure  the  reflected  pressure  pulse  generated  in  the  overburden  water 
during  exploded  cavity  tests.  Measurements  were  made  at  the  bottom  of 
the  pressure  vessel  directly  below  the  charge  as  shown  in  Figure  2.2. 
Although  the  nominal  distance  from  charge  center  to  gage  was  6-5/8  inches 
(17.0  cm),  a  tolerance  of  1/8  inch  (0.318  cm)  should  be  allowed. 

A  summary  of  the  charge  reproducibility  results  is  given  in  Table  2.1, 
where  tests  are  categorized  according  to  the  material  cast  and  the  embedded 
features.  Maximum  reflected  pressure  and  total  reflected  impulse  are 
listed  for  each  test. 

Typical  reflected  pressure  pulses  are  shown  in  Figure  2.12.  The 
corresponding  impulse  curves  are  shown  in  Figure  2.13.  The  pulses  for 
tests  203  and  204  are  associated  with  charge  detonation  in  rock-matching 
grout  (RMG  2C4)  spheres  with  no  embedded  features.  The  average  maximum 
reflected  pressure  for  tests  in  this  category  is  834  psi  (5.75  MPa).  The 
corresponding  average  total  reflected  impulse  is  10.1  psi*msec  (0.070 
MPa*msec).  For  calculation  of  total  impulse,  the  pressure  pulse  was 
assumed  to  end  when  the  pressure  dropped  to  the  initial  value  as  shown 
in  Figure  2.12.  The  pulses  for  tests  215  and  216  are  associated  with 
charge  detonation  in  low  density  rock-matching  grout  (LD  2C4)  spheres 
with  no  embedded  features.  The  average  maximum  reflected  pressure  for 
these  tests  is  250  psi  (17.2  MPa),  which  is  70%  lower  than  the  value  for 
the  RMG  2C4.  The  average  total  reflected  impulse  for  the  LD  2C4  is 
2.8  psi ‘msec  (0.019  MPa*msec),  which  is  73%  lower  than  the  value  for 
RMG  2C4.  Since  the  materials  differ  solely  in  the  air  void  content, 
crushing  of  the  air  voids  in  LD  2C4  is  responsible  for  the  reduction  in 
pressure  and  impulse. 

Table  2.1  shows  that  the  presence  of  a  nearby  vertical  fault  resulted 
in  an  average  maximum  reflected  pressure  of  1088  psi  (7.60  MPa),  which  is 
30%  higher  than  the  value  found  in  tests  without  a  fault.  Similarly,  the 
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average  total  reflected  impulse  for  the  fault  plane  tests  is  13.5  psi*msec 
(0.093  MPa*msec),  an  increase  of  33%.  These  increases  are  attributed  to 
reflection  of  the  pulse  from  the  fault. 

In  conclusion.  Figure  2.12  shows  pulse  shapes  that  are  strikingly 
similar  for  a  given  material,  indicating  a  reproducible  release  of 
energy  from  the  charges. 
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SECTION  3 

HYDROFRACTURE  RESULTS 


3.1  HYDROFRACTURE  TEST  SERIES 

The  experiments  were  performed  to  obtain  hydrofracture  pressures  in 
unexploded  and  exploded  cavities  so  that  by  comparison  the  contribution 
of  the  residual  stress  field  to  containment  may  be  assessed.  They  were 
also  performed  to  provide  the  influence  of  various  factors  on  containment 
such  as  viscosity,  overburden,  rock  properties,  charge  coupling,  and 
nearby  faults  and  tunnels.  Conclusions  based  on  these  results  are 
contained  in  the  summary  of  Section  1.2. 

Results  of  the  hydrofracture  experiments  are  summarized  in  Table  3.1. 
Unexploded  cavity  tests  are  grouped  in  series  1  through  6;  exploded  cavity 

tests  are  grouped  in  series  7  through  17.  Reproducibility  tests,  series  1 

and  7,  were  performed  to  demonstrate  the  reliability  of  the  current 
experimental  techniques  and  to  provide  comparison  data  for  the  following 
areas  of  investigation: 

•  Viscosity  (series  2  and  12).  Dyed  conductive  water, 

with  a  viscosity  of  1  centipoise,  was  chosen  as  the 
standard  hydrofracture  fluid.  To  assess  the  effects 
on  hydrofracture  of  a  more  viscous  fluid,  we  chose 

dyed  glycerol,  with  a  viscosity  of  660  centipoise, 

as  a  substitute  for  water  in  several  unexploded  and 
exploded  cavity  tests. 

•  Overburden  (water)  (series  3).  In  unexploded  cavity  tests, 
overburden  pressure  provides  a  strong  parameter  for  the 
study  of  fracture  initiation  and  growth.  Simulation  of 

of  the  in-situ  stress  state  surrounding  a  typical  under¬ 
ground  nuclear  cavity  was  achieved  by  applying  a  pressure 
of  1000  psi  (6.895  MPa).  Pressures  of  0,  1000,  and  2000 
psi  (0,  6.895,  and  13.79  MPa)  were  found  to  he  sufficient 
for  a  study  of  overburden. 
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•  Overburden  (glycerol)  (series  4).  The  viscosity  and 
overburden  parameters  were  combined  in  several  unexploded 
tests.  Glycerol  hydrofractures  were  performed  on  spheres 
subjected  to  overburden  pressures  of  0,  1000,  and  2000  psi 
(0,  6.895,  and  13.79  MPa). 

•  Material  property  (series  5  and  13) .  Rock-matching  grout 
(RMG.2C4),  which  approximates  the  properties  of  Nevada 
Test  Site  tuff,  was  chosen  as  the  standard  material  tested. 
Effects  of  material  property  variations  on  hydrofracture 
were  studied  by  testing  a  slightly  weaker  low-density  grout 
containing  air  voids  (LD  2C4)  and  a  stronger  granite  simulant 
(GS3) .  Unexploded  and  exploded  cavity  spheres  were  tested. 

Both  water  and  glycerol  hydrofractures  were  performed. 

•  Cavity  lining  (series  6).  A  flexible  cavity  lining  allows 
for  fracture  initiation  and  propagation  without  hydro¬ 
fracture  fluid  acting  directly  on  the  wall  of  the  cavity 
or  entering  the  plane  of  the  crack.  Changes  in  the 
pressures  required  for  fracture  initiation  and  the  subsequent 
stability  of  crack  growth  are  the  major  influences  of  cavity 
lining  on  hydrofracture.  Testing  was  limited  to  unexploded 
cavities. 

•  Vented  cavity  (series  8).  Release  of  the  detonation  products 
from  an  exploded  cavity  before  hydrofracture  provides  the 
limiting  case  of  complete  unloading  for  the  unvented  hydro¬ 
fracture  tests.  Relaxation  of  the  dynamically  produced 
residual  stress  field  surrounding  an  exploded  cavity  is 
enhanced  by  cavity  venting. 

•  Charge  holder  thickness  (series  9).  For  increased  reliability, 
the  charge  holder  wall  thickness  was  increased  from  16  mils  to 
26  mils  (0.406  mm  to  0.660  mm).  All  charges  contained  approxi¬ 
mately  3/8  gram  of  explosive  pressed  to  a  density  of  1  g/eni’. 
The  thinner-walled  charge  produced  a  0.750-inch  (1.905-em) 
exploded  cavity  diameter.  Use  of  the  thicker  walled  charge 
resulted  in  a  0. 800-inch  (2.032-cm)  exploded  cavity  diameter. 

•  Uncoupled  charge  (series  10).  In  standard  exploded  cavity 
tests,  the  explosive  source  is  in  direct  contact  witli  the 
surrounding  material.  In  an  uncoupled  charge  test,  an 
annular  region  of  air  separates  the  explosive  source  from 
the  surrounding  material. 

•  Geometry  (series  11).  In  preparation  for  future  asymmetry 
tests,  several  exploded  cavity  tests  were  performed  on  grout 
cylinders  with  diameters  and  lengths  equal  to  the  diameter 
of  a  sphere  (-12  inches). 
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•  Material  interface  (series  14).  An  important  geological 
feature  at  the  test  site  is  the  presence  of  layers  of 
different  media  in  the  region  surrounding  a  cavity.  This 
feature  was  simulated  in  several  exploded  cavity  hydro¬ 
fracture  tests  by  bonding  rock-matching  grout  to  either  low 
density  rock-matching  grout  or  granite  simulant  along  a 
plane  interface.  The  charge  was  always  detonated  in  the 
rock-matching  grout. 

•  Radial  tunnel  (series  15).  Spheres  containing  hollow  and 
filled  radial  tunnels  were  tested  to  provide  results  for 
possible  tunnel  configurations. 

•  Nearby  tunnel  (series  16).  Spheres  containing  hollow  tunnels 
embedded  at  two  different  distances  from  the  explosive  charge 
were  tested.  The  tunnels  allowed  for  possible  venting  of 
detonation  products  and  loss  of  hydrof racture  fluid. 

•  Fault  plane  (series  17).  Natural  fault  planes  in  the'  vicinity 
of  a  nuclear  device  cavity  were  simulated  by  embedding  sheets 
of  Mylar  in  exploded  cavity  grout  spheres.  Nearby  faults  as 
well  as  a  fault  through  the  explosive  source  were  simulated. 

1.2  UN EXPLODED  CAVTTY  TESTS 
Series  1  -  Reproducibility 

Hydrof racture  tests  were  performed  on  unexploded  cavity  rock-matching 
grout  (RMC  2C4)  spheres  to  demonstrate  the  reliability  of  the  current 
experimental  techniques  and  to  provide  comparison  data  for  parameter 
studies.  The  smooth  unlined  cavity  configuration  shown  in  Figure  2.9  was 
used  in  these  tests.  A  cavity  diameter  of  3/4  inch  (1.90  cm)  was  chosen 
because  it  corresponds  to  the  exploded  cavity  diameter  generated  hv  Lhe 
thinner-walled  3/8-gram  charge.  The  externa]  pressure  applied  to  each 
sphere  was  fixed  at  1000  psi  (6.895  MPa)  to  simulate  the  hydrostatic- 
pressure  in  nuclear  tests  typically  conducted  at  a  depth  of  1100  feet 

3 

(335  m)  in  tuff  with  a  density  of  2.1  g/cm  .  Dyed  conductive  water  was 

3 

pumped  into  each  cavity  at  a  rate  of  4.26  cm  /min.  The  above  parameters 
were  arbitrarily  chosen  as  standard.  The  sensitivity  of  hydrof racture 
pressures  to  variations  in  these  parameters  was  investigated  in  the 
following  series  of  tests. 

Hydrof racture  records  for  tests  93,  94,  160,  and  161  are  shown  in 
Figure  3.1.  Tests  93  and  94  arc  previously  reported  tests  in  which  the 
access  tube  was  not  modified  bv  the  addition  of  the  copper  cone  shown  in 
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Fluid  Detected  at  Surface 


FIGURE  3.1  HYDROFRACTURE  PRESSURES  FOR  UNEXPLODED  CAVITY  TESTS  93,  94,  160, 
AND  161  —  REPRODUCIBILITY 


Figure  2.9.  The  similarity  of  the  four  pressure  records  suggests  that 
this  modification  has  a  negligible  effect  on  hydrofracture  pressures. 

The  results  shown  in  Figure  3.1  are  representative  of  unexploded  cavity 
tests  in  that  a  smooth  increase  in  pressure  is  followed  hv  a  sharp  drop 
to  form  a  well-defined  peak.  Previous  fracture  initiation  studies3  have 
shown  that  the  pronounced  spike  represents  the  initiation  of  a  crack  in 
the  wall  of  the  cavity.  Calibration  tests,  in  which  the  cavity  was 
replaced  by  a  rigid  vessel  of  equivalent  volume,  resulted  in  pressure 
records  witli  an  initial  slope  comparable  to  those  in  a  hvdrofracture  test. 
Hence  tlie  fluid  volume  apparently  required  to  initiate  fracture  can  be 
attributed  to  a  combination  of  fluid  compressibility  and  system  compliance. 
Loss  of  fluid  due  to  porous  flow  is  negligible. 

Fracture  initiation  pressures,  which  provide  the  primary  measure  of 
reproducibility,  were  1950,  1940,  1820,  and  1770  psi  (13.44,  13.38,  12.55, 
and  12.20  MPa).  Reproducibility  is  seen  to  be  excellent,  with  a  maximum 
deviation  of  only  4.3%  from  the  1870-psi  (12.89-MPa)  average. 

Following  fracture  initiation,  a  crack  continues  to  grow  as  cavity 

pressure  decays  slowly  to  a  plateau  that  is  200  to  400  psi  (1.38  to  2.7b 

MPa)  above  the  overburden  pressure.  As  shown  in  Figure  3.1,  hydrofracture 

fluid  advancing  behind  the  crack  tip  was  detected  at  the  surface  of  the 
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sphere  in  test  160  after  6.8  cm  of  flow.  This  volume,  measured  t rom  the 
point  of  fracture  initiation,  provides  an  upper  bound  estimate  of  the 
flow  required  to  generate  a  surface  fracture.  The  subsequent  steady- 
state  cavity  pressure  provides  a  measure  of  the  flow  resistance  along 
the  fracture  plane. 

Figure  3.2  shows  a  typical  fracture  surface  associated  with  an 
unexploded  cavity  test.  The  uniformity  of  the  dyed  region  suggests  that 
the  fracture  initially  propagated  symmetrically  from  the  cavity.  Section¬ 
ing  of  the  unexploded  cavity  spheres  showed  that  a  single  planar  fracture 
surface  was  produced  by  hydrofracLure  in  each  test.  Figure  1.2  also  shows 
a  section  taken  through  the  access  tube  in  test  160.  The  presence  of  the 
copper  cone  prevented  flow  along  the  tube  and  resulted  in  a  horizontal 
fracture  plane.  In  general,  the  cone  allows  for  a  random  orientation  of 
thi'  fracture  plane  in  unexploded  cavity  tests. 
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The  initially  erratic  rise  in  cavity  pressure  in  test  16]  is  attributed 
to  the  presence  of  an  air  void  adjacent  to  the  cavity.  Although  hydro¬ 
fracture  pressures  were  not  significantly  affected  by  the  imperfection, 
the  fracture  plane  did  intersect  this  region.  Previous  tests  with  notched 
cavities^  also  showed  that  imperfections  control  fracture  orientation. 

Series  2  -  Viscosity 

The  effects  of  a  change  in  fluid  viscosity  on  the  pressure  records 
are  shown  in  Figure  3.3.  Tests  160  and  161  are  typical  water  hydro¬ 
fracture  results  selected  from  the  reproducibility  series.  Viscosity 
of  the  dyed  water  was  1  centipoise.  Tile  spheres  in  tests  154  and  155 
are  representative  of  those  hydrof ractured  witii  dyed  glycerol,  which  has 
a  viscosity  of  660  centipoise.  Use  of  the  more  viscous  fluid  increased 
the  average  fracture  initiation  pressure  from  1870  psi  to  2680  psi 
(12.89  MPa  to  18.48  MPa).  For  a  perfectly  smooth  cavity  and  a  completely 
impermeable  sphere,  changes  in  fluid  viscosity  would  have  no  effect  on 
fracture  initiation.  However,  as  fracture  commences  and  fluid  flows 
into  tiie  crack,  the  higher  viscosity  of  the  glycerol  causes  a  spatial 
pressure  decay  from  the  cavity  to  the  crack  tip.  Thus  higher  cavity 
pressures  are  required  to  continue  the  crack  propagation.  A  theoretical 
treatment  of  this  process  is  detailed  in  Appendix  A. 

Figure  3.3  shows  that  the  volume  of  fluid  required  to  establish  a 
steady-state  flow  along  a  fracture  plane  increases  witii  fluid  viscosity. 
However,  the  lower  steady-state  pressures  for  the  glycerol  tests  suggest 
less  resistance  to  flow.  Since  flow  resistance  along  the  planar  crack  is 
controlled  by  crack  width,  use  of  a  more  viscous  hydrofracture  fluid 
apparently  produces  larger  separation  of  the  fractured  surfaces. 

Series  3  -  Overburden  (Water) 

The  effects  of  a  change  in  overburden  pressure  on  water  hydrofracture 
are  shown  in  Figure  3.4.  Tests  165  and  166,  performed  without  external 
pressure,  resulted  in  an  average  fracture  initiation  pressure  of  860  psi 
(5.929  MPa).  Tests  160  and  161,  which  are  typical  results  of  1000-psi 
(6.895  MPa)  overburden,  were  selected  from  the  reproducibility  series. 


Detected  at  Surface 


FIGURE  3.3  HYDROFRACTURE  PRESSURES  FOR  UNEXPLODED  CAVITY  TESTS  154,  155, 
160,  AND  161  —  VISCOSITY 
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Average  fracture  initiation  pressure  for  all  these  tests  was  J870  psi 
(12.89  MPa).  Tests  171  and  172,  performed  at  2000-psi  (13.79-MPa) 
overburden,  resulted  in  an  average  fracture  initiation  pressure  of  2870 
psi  (19.79  MPa).  Results  are  summarized  in  Figure  3.5.  For  water  hydro¬ 
fracture  of  RMG  2C4  spheres,  the  slope  of  the  curve  showing  fracture 
initiation  pressure  versus  overburden  is  one.  The  result  is  consistent 
with  the  fracture  mechanics  analysis  presented  in  Appendix  A,  where 
fracture  initiation  pressure  is  related  to  overburden  by  means  of  fracture 
toughness . 

Figure  3.4  shows  that  an  increase  in  overburden  reduces  the  overall 

rate  of  crack  growth.  With  no  external  pressure  applied,  fluid  was 

3 

detected  at  the  surface  of  a  sphere  after  1.6  cm  of  flow.  With  1000-psi 

3 

(6.895-MPa)  external  pressure,  6.8  cm  of  flow  was  required  before  fluid 
was  detected  at  the  surface.  Hence  overburden  serves  to  increase  fracture 
initiation  pressure  and  slow  crack  growth. 

Series  4  -  Overburden  (Glycerol) 

A  series  of  overburden  tests  was  performed  with  glycerol  as  the 
hydrofracture  fluid  to  provide  additional  data  for  the  study  of  fracture 
initiation  and  growth.  The  pressure  records  are  shown  in  Figure  3.6. 

Test  195,  conducted  with  no  overburden,  resulted  in  a  fracture  initiation 
pressure  of  1760  psi  (12.13  MPa).  Average  fracture  initiation  pressure 
for  tests  154  and  155,  conducted  at  1000-psi  (6.895-MPa)  overburden,  was 
2580  psi  (17.79  MPa).  Finally,  test  178,  conducted  with  2000-psi  (13.79-MPa) 
overburden,  resulted  in  cavity  fracture  at  3610  psi  (24.89  MPa).  A  summary 
of  these  results  is  shown  in  Figure  3.5.  For  glycerol  hydrofracture  of 
RMG  2C4  spheres,  the  slope  of  the  curve  showing  fracture  initiation 
pressure  versus  overburden  is  one.  The  result,  which  was  also  obtained 
for  water  hydrofracture,  is  consistent  witli  the  fracture  mechanics  analysis 
presented  in  Appendix  A. 

The  effect  of  hydrofracture  fluid  viscosity  on  crack  growth  and 

fluid  motion  is  shown  by  a  comparison  of  the  zero  overburden  tests  for 

water  and  glycerol.  In  these  tests  visual  detection  of  surface  cracking 

and  fluid  arrival  was  convenient.  For  glycerol  hydrofracture.  Figure  3.6 
3 

shows  that  2.1  cm  of  flow  was  required  to  produce  surface  cracking, 
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FIGURE  3.5  FRACTURE  INITIATION  PRESSURE  VERSUS  OVERBURDEN  FOR  UNEXPLODED 
CAVITY  TESTS 
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and  b.l  cm  of  flow  was  required  to  drive  fluid  to  the  surface, 
measurements  are  made  relative  to  the  point  of  fracture  initiation. 

For  water  hydrofracture.  Figure  3.4  shows  that  an  average  flow  of  1.5  cin^ 
was  required  to  drive  fluid  to  the  surface.  Cracking  of  the  surface 
was  not  detected  before  fluid  arrival  because  of  the  initially  small 
crack  width.  The  conclusion  drawn  from  the  data  is  that  lower  viscosity 
fluid  flows  nearer  the  crack  tip  and  generates  increased  rates  of  crack 
growth  and  fluid  motion. 

Series_ 5  -  Material  Property 

Kffects  of  material  properties  on  containment  were  studied  by 
comparing  hydrofracture  records  for  rock-matching  grout  (RMG  2C4),  low 
density  rock-matching  grout  (LD  2C4),  and  granite  simulant  (GS3)  spheres. 
Properties  of  these  materials  are  given  in  Appendix  C. 

Hydrofrac.ture  pressures  for  LD  2C4  tests  193  and  194  are  shown  in 
Figure  3.7.  The  records  have  overall  features  typical  of  unexploded 
cavity  tests,  but  differ  in  detail.  For  example,  the  rise  in  pressure 
immediately  following  cavity  fracture  indicates  a  temporarily  stable 
crack  growth.  The  feature  is  atypical  of  unexploded  cavity  tests  but 
has  been  observed  in  tests  on  RMG  2C4  and  GS3.  The  average  fracture 
initiation  pressure  is  1580  psi  (10.89  MPa).  Hydrofracture  records  for 
tests  lbO  and  161  were  selected  from  the  reproducibility  series  and  are 
also  shown  in  Figure  3.7  for  comparison.  The  average  fracture  initiation 
pressure  for  RMG  2C4  tests  is  1870  psi  (12.89  MPa).  Hence  the  average 
fracture  initiation  pressure  for  LD  2C4  is  15%  less  than  the  value  for 
RMG  2C4.  A  decrease  is  expected  since  the  tensile  strength  of  LD  2C4 
[460  psi  (3.172  MPa)]  is  13%  less  than  the  tensile  strength  of  RMG  2G4 
[530  psi  (3.654  MPa)],  which  implies  a  lower  fracture  toughness. 

Pressure  records  for  tests  on  GS3  spheres  are  shown  in  Figure  3.8. 
These  results  show  the  effects  of  overburden  and  fluid  viscosity  on 
hydrofracture.  Tests  224  and  222  are  water  hydrofractures  conducted  at 
0-psi  and  1000-psi  (6.895-flPa)  overburden,  respectively.  The  corres¬ 
ponding  fracture  initiation  pressures  are  1940  psi  (13.36  MPa)  and 
2760  psi  (19.03  MPa).  Tests  196  and  218  are  glycerol  hydrof t a  lures 
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FIGURE  3.7  HYDROFRACTURE  PRESSURES  FOR  UNEXPLODED  CAVITY  TESTS  160,  161,  193, 
AND  194  —  MATERIAL  PROPERTY  (LD  2C4) 
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AND  224  —  MATERIAL  PROPERTY  (GS3) 


conducted  at  0-psi  and  1000-psi  (6 . 89  5-MPa )  overburden,  respectively. 

The  corresponding  fracture  initiation  pressures-:  are  3231)  psi  (22.27  MPa) 
and  4840  psi  (33.37  MPa).  These  values  are  shown  in  Figure  3.5.  For 
water  hydrofracture,  tie  slope  of  the  curve  showing  fracture  initiation 
versus  overburden  is  one.  This  result,  which  agrees  wit!  the  results  of 
water  and  glycerol  hydrofracture  tests  on  RMC  2C4,  supports  the  conclusion 
of  Appendix  A  regarding  the  dependence  of  fracture  initiation  pressure  on 
overburden.  For  glycerol  hydro!" raeture  of  the  two  (IS3  spheres,  the  slope 
of  the  curve  showing  fracture  initiation  pressure  versus  overburden  is 
greater  titan  one.  However,  if  a  typical  scatter  of  unexploded  cavity 
results  is  considered,  the  data  points  can  be  fitted  to  a  line  with  unit 
slope.  Since  the  tensile  strength  of  C1S3  [920  psi  (6.343  MPa)]  is  74/ 
greater  than  the  tensile  strength  of  KMC  2C4  [530  psi  (3.654  MPa)],  implying 
greater  fracture  toughness,  the  increased  fracture  initiation  pressures 
associated  with  CS3  are  expected. 

Series  _6  _-__Cavitj  Lin  wig 

Unexploded  cavities  were  lined  with  an  impermeable  membrane  using 
the  technique  described  in  Section  2.2.  The  lining  is  intended  to 
prevent  hydrofracture  fluid  from  acting  directly  on  the  cavity  wall 
and  from  entering  the  fracture  plane  following  cavity  fracture.  A 
comparison  of  pressure  records  for  lined  and  unlined  cavities  then 
allows  assessment  of  the  effects  of  hydro!" raeture  fluid  on  fracture 
initiation  pressure  and  stability  of  subsequent  crack  growth.  Crack 
growth  is  stable  if  increasing  cavity  pressure  is  required  to  propagate 
the  crack. 

Water  hydrofracture  records  for  lined  cavity  tests  179  and  ISO  are 
shown  in  Figure  3.9.  In  each  test  a  small  leak  path  developed  around 
the  access  tube  before  fracture  initiation  that  allowed  tluid  to  escape. 

Loss  of  fluid  from  the  cavity  lining  is  indicated  by  the  small  pressure 
spikes  on  the  hvdrof rac tur  •  record  below  2000  ps i  (13. 7l)  MPa).  With  the 
restricted  flow,  an  average  fracture  initiation  pressure  of  2520  psi 
(17.37  MPa)  was  developed.  Fracture  initiation  is  indicated  by  the 
pronounced  pressure  spikes.  In  each  test  stable  crack  growth  followed 
fracture  initiation.  In  test  180,  a  cavity  pressure  of  5060  psi  (35.09  MPa) 
was  .applied  without  generating,  a  surface  crack.  Pumping  was  then  stopped. 
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HYDROFRACTURE  PRESSURES  FOR  UNEXPLODED  CAVITY  TESTS  160,  161, 
179,  AND  180  —  CAVITY  LINING 


Subsequent  inspection  of  the  sphere  revealed  that  hydrofracture  fluid 

flowed  5  inches  (12.7  cm)  from  the  center  of  the  sphere.  In  test  179, 

cavity  pressure  reached  3910  psi  (26.96  MPa)  before  crack  growth  became 

3 

unstable.  Surface  cracking  was  detected  after  6.5  cm  of  flow.  Hydro¬ 
fracture  records  for  equivalent  unlined  cavity  tests  160  and  161  are 
shown  in  Figure  3.9  for  comparison.  These  tests  were  selected  from  the 
reproducibility  series  where  the  average  fracture  initiation  pressure 
was  1870  psi  (12.89  MPa).  Hence  the  partial  cavity  lining  increased  the 
average  fracture  initiation  pressure  by  35 7. 

3.3  EXPLODED  CAVITY  TESTS 
Series  7  -  Reproducibility 

Unvented  exploded  cavity  hydrofracture  tests  were  performed  to  assess 
the  reproducibility  of  the  experimental  techniques  and  to  provide  reference 
data  for  the  parameter  studies.  The  configuration  and  procedure  for 
conducting  these  tests  were  described  in  Section  2.2.  With  this  configu¬ 
ration  the  cavity  gas  pressure  resulting  from  charge  detonation  was 
measured.  Also,  pumping  was  started  at  the  time  of  charge  detonation 
so  that  stress  relaxation  due  to  loss  of  cavity  pressure  and  creep  could 
be  minimized.  The  compliance  of  the  system  was  reduced  by  prepressurizing 
the  hydrofracture  fluid  to  a  nominal  value  of  300  psi  (2.068  MPa). 

Hydrofracture  records  for  tests  158,  159,  169,  and  170  are  shown 
in  Figure  3.10.  Not  shown  are  the  initially  high  cavity  gas  pressures 
that  prohibited  the  immediate  lowering  of  the  steel  ball  al  the  end  of 
the  access  tube  (Figure  2.7).  As  the  gas  cooled,  the  pressure  decayed  and 
allowed  the  access  tube  to  be  opened.  The  cavity  gases  then  reached  a 
state  of  equilibrium  with  the  hydrofracture  fluid.  This  process,  which 
required  approximately  four  seconds  to  complete,  precedes  the  liv'd  ro- 
Iracture  records  shown  in  Figure  3.10.  A  theoretical  estimate  ol  the 
gas  pressure  in  a  dynamically  expanding  cavity  is  shown  in  Figure  3.11. 

The  curve  was  derived  from  the  pressure-volume  data  listed  in  Table  3.2. 

The  explosive  source  used  in  the  calculations  was  assumed  to  be  a  mixture 
of  the  PETN  and  Lucite  composing  the  charge  shown  in  Figure  2.5.  The 
assumption  is  based  on  the  experimental  observation  that  the  Lucite 
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FIGURE  3.10  HYDROFRACTURE  PRESSURES  FOR  UNVENTED  EXPLODED  CAVITY  TESTS 
158,  159,  169,  AND  170  —  REPRODUCIBILITY 
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FIGURE  3.11  CALCULATED  CAVITY  GAS  PRESSURE  FOR  A  DYNAMICALLY 
EXPANDING  CAVITY 
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PRESSURE  IN  EXPANDING  CAVITY  x  10'3  —  MPa 
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Pressure-Volume  data  from  SRI  TIGER  code  (PETN  density  1  gram/cm3, 
BKW  equation  of  state). 

Pinal  cavity  radius  ratio  for  GS3  (interpolated). 

Final  cavity  radius  ratio  for  RMG  2CU . 

Final  cavity  radius  ratio  for  U>  2C4  (interpolated). 


charge  holder  is  completely  consumed  in  the  explosion.  Since  the  exploded 
cavity  diameter  is  twice  the  initial  cavity  diameter,  the  theoretical 
estimate  of  cavity  pressure  immediately  following  detonation  is  found  to 
he  10,100  psi  (75.84  MPa).  This  is  the  pressure  that  initially  holds  the 
steel  ball  against  the  end  of  the  access  tube. 

Figure  3.10  shows  that  hydrofracture  pressures  for  unvested  exploded 
cavity  tests  are  related  to  the  unloading  of  cavity  pressure.  In  test 
158,  pressure  decayed  to  1110  psi  (7.653  MPa),  and  the  subsequent  pressure 
required  to  initiate  cavity  fracture  was  3140  psi  (21.65  MPa).  In  test 
170,  a  minimum  cavity  pressure  of  1460  psi  (10.07  MPa)  resulted  in  a 
fracture  initiation  pressure  of  3530  psi  (24.34  MPa).  Despite  these 
differences,  reproducibility  is  good.  The  maximum  deviation  from  the 
3290-psi  (22.68-MPa)  average  fracture  initiation  pressure  is  8.5%.  This 
average  is  76%  higher  than  the  corresponding  1870-psi  (12.89-MPa)  average 
for  equivalent  unexploded  cavity  tests. 

The  existence  of  a  compressive  residual  stress  field  surrounding  an 
exploded  cavity  is  indicated  by  the  high  hvdrofracture  pressures  and  the 
general  stability  of  early  crack  growth  shown  in  Figure  3.10.  in  addition. 
Figure  3.12  shows  an  undyed  region  in  the  neighborhood  of  the  cavity  in 
test  170.  Hydrofracture  fluid  appears  to  have  flowed  from  the  eavitv  along 
isolated  paths  and  spread  around  an  annular  region  of  compressive  residual 
stress . 

Series  8  -  Vented  Cavity 

Vented  exploded  cavity  tests  were  performed  to  assess  the  effect  on 
containment  of  complete  unloading  of  cavity  gas  pressure.  The  configuration 
and  procedure  for  conducting  these  tests  were  described  in  Section  2.2.  In 
each  test,  water  hydrofracture  was  startl’d  20  to  40  minutes  following  charge 
detonation.  The  analysis  presented  in  Appendix  li  indicates  that  relaxation 
of  the  explosively  generated  stress  field  during,  this  time  interval  caused 
hv  pore  water  migration  may  he  negligible.  However,  creep  may  give  rise 
to  different  residual  stress  fields  by  the  time  pumping,  begins. 
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Hydrofracture  records  for  tests  162,  163,  164,  and  177  are  si  town  in 
Figure  3.13.  In  test  164,  cavity  pressure  increased  smoothly  to  a  fracture 
initiation  spike  and  then  decreased  monotonical ly  to  a  steady  state.  In 
the  remaining  tests,  fracture  initiation  was  followed  by  a  small  increase 
in  pressure  before  the  monotonic  decay  to  steady  state.  Stability  of 
crack  growth  thus  appears  to  be  marginal.  However,  compared"  with  the 
initially  rapid  decay  in  cavity  pressure  for  unexploded  cavity  tests,  the 
initially  slow  decay  in  cavity  pressure  for  the  exploded  cavity  tests 
indicates  that  compressive  stresses  are  tending  to  stabilize  crack  growth. 

In  addition,  the  average  steady-state  cavity  pressure  for  the  exploded 
cavity  tests  is  500  psi  (3.45  MPa)  higher  than  the  value  for  the  corres¬ 
ponding  unexploded  cavity  tests.  Since  steady-state  pressure  is  a  measure 
of  flow-resistance  along  a  fracture  plane,  a  higher  value  suggests 
compressive  stresses  restricting  crack  growth.  Figure  3.14,  which  shows 
the  major  fracture  surface  for  test  164,  reveals  that  hydrcf racture  fluid 
escaped  from  the  cavity  along  a  fracture  plane  of  small  sectorial  extent 
and  then  spread  around  the  annular  region  of  residual  compressive  stress. 

Fracture  initiation  pressures  for  the  vented  exploded  cavity  tests 
ranged  from  2080  to  2930  psi  (14.34  to  20.20  MPa).  Values  at  the  lower 
end  of  the  range  are  slightly  higher  than  those  for  the  corresponding 
unexploded  cavity  tests.  Values  at  the  upper  end  of  the  range  are  slightly 
lower  than  the  minimum  fracture  initiation  pressure  for  the  corresponding 
unvented  exploded  cavity  tests. 

Series  9  -  Charge  Holder  Th ickness 

Consistent  detonation  of  the  explosive  source  cast  in  a  sphere  requires 
that  the  Lucite  charge  holder  carry  the  stresses  associated  with  overburden. 
Although  the  design  shown  in  Figure  2.5  is  generally  adequate,  reliability 
was  improved  by  increasing  the  outside  diameter  of  the  spherical  charge 
holder  from  0.375  inch  (0.952  cm)  to  0.400  inch  (1.016  cm).  The  improved 
charge  design,  shown  in  Figure  2.6,  was  incorporated  in  all  exploded 
cavity  tests  numbered  199  and  higher. 
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Figure  3.15  shows  hydrofracture  records  for  unvent ed  exploded  cavity 
tests  199,  200,  and  204  performed  with  the  improved  charge.  Shown  in  the 
figure  for  comparison  are  records  from  the  reproducibility  series  (tests 
159  and  170)  obtained  with  the  thinner-walled  charge.  Hydrofracture  tests 
conducted  with  the  new  design  yielded  pressures  within  the  range  of  pressures 
from  tests  with  the  thinner-walled  charge.  Although  the  new  charge  results 
in  a  slightly  larger  exploded  cavity  diameter  [0.800  inch  (2.032  cm) 
compared  with  0.750  inch  (1.905  cm)],  the  ratio  of  exploded  cavity  diameter 
to  charge  diameter  remains  a  factor  of  two.  Figure  3.16  shows  the  major 
fracture  surface  associated  with  test  200.  The  dyed  pattern  is  typical 
of  unvented  exploded  cavity  tests.  Fluid  escaped  from  the  cavity  along 
an  area  of  small  sectorial  extent  and  then  spread  around  an  annular  region 
of  compressive  residual  stress. 

The  pressure  record  for  test  199  reflects  the  importance  of  adequately 
pressurizing  the  hydrofracture  system  before  charge  detonation.  In  this 
test  the  nominal  300-psi  (2.068-MPa)  prepressure  was  eliminated.  System 
compliance  with  no  prepressure  required  15  seconds  of  pumping  to  open  the 
access  tube.  For  a  prepressure  of  300  psi  (2.068  MPa),  only  3  to  5  seconds 
of  pumping  is  normally  required.  Hence  the  hydrofracture  record  for  test 
199  begins  15  seconds  after  detonation  with  an  initial  cavity  pressure  of 
970  psi  (6.688  MPa).  The  subsequent  fracture  initiation  pressure,  3080  psi 
(21.24  MPa),  is  at  the  low  end  of  the  range  for  unvented  exploded  cavity 
tests . 


Series  10  -  Uncoupled  Charge 

Water  hydrofracture  tests  were  performed  on  unvented  exploded  cavity 
spheres  in  which  a  3/16-inch  (0.476-cm)  annular  air  space  separated  the 
charge  from  the  surrounding  grout  before  detonation.  This  configuration, 
shown  schematically  in  Figure  3.17,  simulates  an  uncoupled  nuclear  cavity 
detonation.  The  initial  3/4-inch-diameter  (1.90 -cm)  cavity  was  cast 
around  a  spherical  glass  shell  with  an  8-mil  (0.203-mm)  wall  thickness. 
The  shell  shattered  upon  charge  detonation. 
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Hydrofracture  records  for  uncoupled  charge  tests  174,  175,  and  176 
are  shown  In  Figure  3.17.  For  comparison,  tests  159  and  170  from  the 
reproducibility  series  are  also  shown.  Uncoupling  the  charge  allows  the 
detonation  products  to  expand  into  a  larger  volume.  Hence  unloading  of 
cavity  gas  pressure  is  more  pronounced  in  the  uncoupled  charge  configuration. 
The  exploded  cavity  diameters  in  tests  174,  175,  and  176  were  54/64,  61/64, 
and  62/64  inch  (2.14,  2.42,  and  2.46  cm),  respectively.  Expansion  of 
cavity  diameters  thus  ranged  from  a  minimum  of  6/64  to  a  maximum  of 
14/64  inch  (0.238  to  0.556  cm).  The  highest  residual  cavity  gas  pressure, 

1020  psi  (7.033  MPa),  and  the  highest  fracture  initiation  pressure, 

3720  psi  (26.65  MPa),  for  the  uncoupled  charge  experiments  were  developed 
during  test  174,  where  cavity  expansion  was  the  smallest. 

In  tests  175  and  176,  embedded  fluid  arrival  gages  detected  hydro¬ 
fracture  fluid  1  inch  (2.54  cm)  from  the  center  of  the  sphere  after 
3  3 

0.92  cm  and  2.18  cm  of  flow,  respectively.  Since  cavity  pressure  was 
increasing  at  this  time,  crack  growth  was  stable.  In  test  174,  fluid  was 

3 

detected  2  inches  (5.08  cm)  from  the  center  of  the  sphere  after  0.97  cm 
of  flow.  Since  cavity  pressure  was  unchanging  at  this  time,  crack  growth 
was  neutral. 

Figure  3.18  shows  the  multiple  fracture  planes  that  developed  during 
hydrofracture  in  test  176.  A  complex  pattern  of  fracture  surfaces  is  a 
general  characteristic  of  exploded  cavity  tests. 

Series  11  -  Geometry 

To  study  the  influence  of  different  external  geometries  on  hydro¬ 
fracture  pressures  and  for  use  in  future  asymmetry  tests,  we  performed 
water  hydrofracture  tests  181  and  132  on  cylinders  with  unvented  exploded 
cavities.  The  cylinders  measured  12  inches  (30.48  cm)  in  diameter  and 
12  inches  (30.48  cm)  in  length.  Hydrofracture  records  are  shown  in 
Figure  3.19,  where  a  comparison  is  made  with  typical  tests  conducted  on 
12-inch-diameter  (30.48-cm)  spheres.  Although  the  records  for  the  cylinder 
exhibited  no  new  features,  the  change  in  geometry  reoriented  the  major 
fracture  plane.  For  a  sphere,  surface  cracking  generally  appears  in  a 
plane  inclined  to  the  access  tube.  For  a  cylinder,  a  vertical  fracture 
plane  was  developed  as  shown  in  Figure  3.20. 
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FIGURE  3.19  HYDROFRACTURE  PRESSURES  FOR  UNVENTED  EXPLODED  CAVITY  TESTS 
159,  170,  181,  AND  182  —  GEOMETRY 


Scries  12-  Viscosity 

Glycerol  hydrofracture  of  unvented  exploded  cavity  spheres  was 
performed  to  study  the  influence  of  fluid  viscosity  on  hydrofracture. 
Pressure  records  for  tests  135,  142,  and  157  are  shown  in  Figure  3.21. 

For  comparison,  typical  water  hydrofracture  tests  159  and  170  are  also 
shown.  Tests  135  and  142  are  previously  reported  experiments ^  in  which 
cavity  gas  pressure  was  allowed  to  decay  to  near  steady  state  before 
fluid  was  pumped  into  the  cavity.  Test  157  was  performed  during  the 
developmental  stage  of  unvented  cavity  tenting  when  the  delay  from  charge 
detonation  to  measurement  of  cavity  gas  pressure  was  longer  than  the 
nominal  4-second  delay  in  later  tests.  Hence  the  full  effect  of  a  change 
in  fluid  viscosity  on  hydrofracture  is  not  indicated  in  Figure  3.21.  A 
comparison  of  the  hydrofracture  pressures  for  tests  157  and  159  does 
indicate  the  trend,  however.  Unloading  of  cavity  gas  pressure  was  the  same 
in  both  tests.  For  water  hydrofracture,  the  fracture  initiation  pressure 
was  3010  psi  (20.75  MPa).  For  glycerol  hydrofracture,  the  fracture 
initiation  pressure  was  3540  psi  (24.41  MPa),  an  increase  of  18%. 

Series  13  -  Material  Property 

Effects  of  a  change  in  material  properties  on  containment  were 
evaluated  by  comparing  unvented  exploded  cavity  hydrofracture  records 
for  rock-matching  grout  (RMG  2C4) ,  low  density  rock-matching  grout  (LD  2C4), 
and  granite  simulant  (GS3)  spheres.  Properties  of  these  materials  are 
presented  in  Appendix  C. 

Figure  3.22  shows  the  water  hydrofracture  records  for  unvented 
exploded  cavity  tests  189,  215,  and  216  performed  on  LD  2C4  spheres.  For 
comparison,  results  of  tests  159  and  170  performed  on  RMG  2C4  spheres  are 
also  shown.  An  important  feature  of  the  low  density  material  in  exploded 
cavity  tests  is  the  13.4%  air  void  content.  This  property  accounts  for 
the  formation  of  an  exploded  cavity  with  a  final-to-in i t ial  radius  ratio 
of  2.4.  Figure  3.11  shows  that  the  corresponding  theoretical  estimate  of 
initial  gas  pressure  in  the  expanded  cavity  is  4300  psi  (29.65  MPa).  For 
a  RMG  2C4  sphere,  the  radius  ratio  is  2.0  and  the  theoretical  initial 
pressure  is  10,100  psi  (69.64  MPa).  As  shown  in  Figure  3.22,  these  values 
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are  reflected  In  the  initial  cavity  pressures  and  the  volumes  of  fluid 
required  to  fill  the  cavities.  The  average  fracture  initiation  pressure 
for  the  l.D  2C4  snheres,  1770  psi  (12.20  MPa),  is  considerably  less  than 
the  3290-psi  (22.68-MPa)  average  for  RMG  2C4  spheres,  but  is  comparable 
to  the  1580-psi  (10.89-MPa)  average  for  unexploded  cavity  low  density 
spheres.  Hence,  the  residual  stress  field  developed  around  an  LD  2C4 
exploded  cavity  is  weak.  Figure  3.23  shows  the  fracture  surface  for  test 
216.  The  absence  of  a  clearly  defined  undyed  region  surrounding  the  cavity 
is  further  evidence  of  weak  residual  stresses. 

Dynamic  pressure  pulses  measured  at  the  bottom  of  the  containment 
vessel  were  shown  in  Figure  2.12.  For  LD  2C4  the  average  reflected  peak 
pressure  is  250  psi  (1.72  MPa),  and  the  average  total  reflected  impulse 
is  2.8  psi*msec  (0.019  MPa*msec).  The  corresponding  values  for  RMG  2C4 
tests  are  834  psi  (5.75  MPa)  and  10.1  psi*msec  (0.070  MPa*msec).  Hence 
attenuation  of  the  stress  wave  associated  with  charge  detonation  is  more 
pronounced  in  the  low  density  material.  Explosive  crushing  of  the  air 
voids  contained  in  LD  2C4  is  the  process  by  which  energy  is  absorbed  without 
the  generation  of  significant  residual  stresses. 

The  results  of  additional  tests  on  materials  containing  air  voids 
are  shown  in  Figure  3.24.  Tests  201  and  202  were  unvented  exploded  cavity 
water  hydrofractures  performed  on  a  medium  density  rock-matching  grout 
(MD  2C4)  containing  approximately  6%  air  voids.  For  comparison,  results 
for  the  fully  saturated  RMG  2C4  and  the  LD  2C4  witli  13.4%  air  voids  are 
also  shown.  Tiie  values  for  exploded  cavity  diameter  and  fracture  initiation 
pressure  obtained  in  the  MD  2C4  tests  lie  between  the  corresponding  values 
obtained  in  the  RMG  2C4  and  LD  2C4  tests.  Hence  the  air  void  content  of  a 
material  represents  a  significant  containment  parameter. 

A  final  series  of  material  property  hydrofracture  tests  was  performed 
on  granite  simulant  (GS3)  spheres.  Charge  detonation  resulted  in  a  final- 
to-initial  cavity  radius  ratio  of  1.3.  Figure  3.11  shows  that  the  corres¬ 
ponding  theoretical  estimate  of  initial  cavity  pressure1  is  97,400  psi 
(67.15  MPa).  Tn  the  tests,  high  gas  pressures  resulted  in  fracture 
initiation  before  fluid  was  pumped  into  the  cavity. 
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FIGURE  3.24  HYDROFRACTURE  PRESSURES  FOR  UNVENTED  EXPLODEC 
170,  201,  202,  AND  216  —  MATERIAL  PROPERTY  (MD  2C4) 


Series  14  -  Material  Interface 

Potential  sites  for  underground  nuclear  testing  include  areas  near 
a  surface  of  discontinuity  that  separates  regions  of  different  materials. 

The  effect  of  this  geological  feature  on  containment  was  simulated  by 
conducting  unvented  exploded  cavity  tests  on  spheres  containing  a 
material  interface. 

The  case  of  a  weaker  nearby  material  was  examined  by  testing  composite 
spheres  of  rock-matching  grout  (RMG  2C4)  and  low  density  rock-matching 
grout  (LD  2C4).  In  each  test  the  plane  horizontal  interface  passed  within 
3/4  inch  (1.905  cm),  or  approximately  one  exploded  cavity  diameter,  of 
the  center  of  the  sphere.  The  locations  of  the  interface  and  the  resulting 
hydrofracture  records  for  tests  205  and  208  are  shown  in  Figure  3.25.  Charge 
detonation  occurred  in  the  RMG  2C4,  and  the  resulting  exploded  cavity 
diameter,  0.800  inch  (2.032  cm),  was  the  same  as  that  found  in  tests  on 
VMG  2C4  spheres.  Hydrofracture  produced  a  fracture  plane  that  propagated 
across  the  interface  without  changing  direction.  Figure  3.25  shows  that 
hydrol racture  pressures  for  the  composite  spheres  are  similar  to  those 
for  standard  RMG  2C4  tests.  Hence  the  influence  of  the  nearby  low-density 
material  appears  to  be  negligible. 

The  case  of  a  stronger  nearby  material  was  examined  by  testing  composite 
spheres  of  RMG  2C4  and  granite  simulant  (GS3) .  The  distance  of  the  inter¬ 
face  from  the  center  of  the  sphere  was  2  inches  (5.080  cm)  in  test  188  and 
3/4  inch  (1.905  cm)  in  test  213,  as  shown  in  Figure  3.26.  Charge  detonation 
again  occurred  in  RMG  2C4  and  resulted  in  an  exploded  cavity  diameter  the 
same  as  that  found  in  tests  on  2C4  spheres.  The  pressure  record  for  test 
188  is  typical  of  a  test  on  a  2C4  sphere,  with  fracture  initiation  at 
3540  psi  (24.41  MPa).  The  pressure  record  for  test  213,  however,  indicates 
fracture  initiation  at  1750  psi  (12.7  MPa).  The  low  fracture  initiation 
pressure  suggests  low  adherence  of  RMG  2C4  to  GS3.  In  fact,  the  hydrofractur 
driven  crack  propagated  along  the  interface.  Bonding  in  test  188  appeared 
to  be  adequate,  however.  Hence  boundary  conditions  along  an  interface 
appear  to  control  hydroi racture  pressures  in  the  composite  spheres  tested. 
Future  composite  sphere  experiments  with  GS3  requirt'  development  of  a 
technique  for  producing  a  reliable  interface  bond. 
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FIGURE  3.25  HYDROFRACTURE  PRESSURES  FOR  UNVENTED  EXPLODED  CAVITY  ' 
159.  170,  205,  AND  208  —  MATERIAL  INTERFACE  (RMG  2C4/LD  2C4) 
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IGURE  3.26  HYDROFRACTURE  PRESSURES  FOR  UNVENTED  EXPLODED  CAVITY  TESTS 
159.  170,  188,  AND  213  —  MATERIAL  INTERFACE  (RMG  2C4/GS3) 


Series  15  -  Radial  Tunnel 


An  underground  nuclear  test  has  a  tunnel  extending  from  the  device 
cavity.  The  effects  of  tunnels  on  containment  were  studied  by  conducting 
unvented  exploded  cavity  tests  on  spheres  containing  different  radial 
tunnel  configurations. 

Simulation  of  a  radial  tunnel  filled  over  most  of  its  length  was 

provided  by  the  configuration  shown  in  Figure  3.27(a).  A  solid  hucite 

rod  filled  the  entire  length  of  tunnel  except  for  a  1/8-inch  (0.318-cm) 

air  gap  at  the  end  opposite  the  charge.  The  resulting  hydrofracture 

record  for  test  185  is  shown  in  Figure  3.28.  The  cavity  pressure  and 

maximum  hydrofracture  pressure  were  consistent  with  those  in  tests  159 

and  170  without  a  tunnel.  Examination  of  the  sphere  following  the  test 

revealed  that  the  Lucite  rod  had  been  extruded  to  the  end  of  the  tunnel. 

However,  the  fracture  plane  reaching  the  surface  did  not  intersect  the 

3 

tunnel.  Surface  cracking  was  detected  after  3.8  cm  of  flow,  with 
volume  measured  from  the  point  of  fracture  initiation.  The  associated 
sudden  loss  of  cavity  pressure  shown  in  Figure  3.28  indicates  development 
of  widespread  surface  cracking. 

Simulation  of  a  hollow  radial  tunnel  was  provided  by  the  configura- 
tionshownin  Figure  3.27(b).  The  2-1  /2-inch-long  (6.35-cin)  tunnel  extend¬ 
ing  from  the  charge  was  cast  around  a  thin-walled  glass  tube,  which 
shattered  upon  charge  detonation.  Hence,  the  pressure  record  for  test 
214,  shown  in  Figure  3.28,  is  the  hydrofracture  associated  with  a  hollow 
unlined  tunnel.  For  comparison,  results  of  tests  without  a  tunnel,  159 
and  170,  are  also  shown.  Since  venting  of  the  cavity  gases  to  the  open 
tunnel  reduced  the  initial  cavitv  gas  pressure,  the  subsequent  hydro¬ 
fracture  pressures  were  also  reduced.  Fracture  initiated  from  the  end  of 
the  tunnel  when  the  cavity  pressure  reached  1560  psi  (10.76  MPa).  Examina¬ 
tion  of  the  sphere  following  the  test  showed  no  signs  of  stemming  along  the 
tunne 1 . 
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FIGURE  3.27  CONFIGURATIONS  FOR  RADIAL  TUNNEL  TESTS 
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FIGURE  3.28  HYDROFRACTURE  PRESSURES  FOR  UNVENTED  EXPLODED  CAVITY  TESTS 
159,  170,  185,  AND  214  —  RADIAL  TUNNEL 


Series  16  -  Nearby  Tunnel 

Nearby  tunnels  associated  with  past  nuclear  tests  represent  a 
possible  path  for  the  venting  of  detonation  products  in  future  tests. 

The  effects  of  these  test  site  features  on  containment  were  examined  by 
testing  unvented  exploded  cavity  spheres  containing  nearby  tunnels. 

The  configuration  for  simulating  a  nearby  tunnel  is  shown  in 
Figure  3.29.  A  thin-walled  plastic  tube  serves  as  a  lining  for  the 
6-inch-long  (15.2-cm)  tunnel.  The  distance  d  from  the  center  of  the 
tunnel  to  the  center  of  the  charge  provides  a  parameter  for  locating  the 
tunnel . 

In  test  186  the  distance  d  was  1/2  inch  (1.27  cm).  Venting  of  the 
explosive  gases  into  the  tunnel  upon  charge  detonation  resulted  in  the 
low  residual  cavity  gas  pressure  measurement  shown  in  Figure  3.30. 

Subsequent  hydrofracture  pressures  were  also  low  relative  to  those  for 
tests  159  and  170  without  a  tunnel.  Hydrofracture  resulted  in  fluid 
flowing  from  the  cavity  into  the  tunnel  along  the  explosively  formed 
crack.  Fluid  then  flowed  along  the  tunnel  and  was  discharged  along  a 
1-1/2-inch-long  (3.81-cm)  section  at  one  end  of  the  tunnel.  As  shown  in 
Figure  3.31,  a  3-inch-diameter  (7.62-cm)  region  surrounding  the  cavity 
was  left  undyed.  This  undyed  region  is  further  evidence  of  the  residual 
stress  field  surrounding  an  exploded  cavity.  The  absence  of  a  pronounced 
fracture  initiation  spike  in  the  pressure  record  for  test  186  confirms  that 
the  cavity  ruptured  as  a  result  of  charge  detonation. 

In  tests  209  and  223,  the  distance  d  was  3/4  inch  (1.90  cm),  or 
approximately  one  exploded  cavity  diameter.  The  detonation  products  were 
contained  in  the  cavity  and  generated  iiydrof racture  records  typical  of  tests 
without  a  tunnel  (Figure  3.30).  In  both  tests  the  crack  extending  to  the 
surface  of  the  sphere  intersected  the  tunnel  but  did  not  rupture  the 
plastic  liner.  Hence  the  tunnels  appeared  to  have  no  influence  on  hydro¬ 
fracture.  The  sudden  drop  in  cavity  pressure  shown  in  Figure  3.30  for 
test  223  corresponds  to  the  detection  of  a  surface  crack. 
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FIGURE  3.29  CONFIGURATION  FOR  NEARBY  TUNNEL  TESTS 
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Series  17  -  Fault  Plane 


Fault  planes  in  tlie  vicinity  of  a  nuclear  device  cavity  represent 
paths  along  which  the  detonation  products  may  escape  to  the  atmosphere. 
The  effects  of  these  geological  features  on  containment  were  examined  by 
testing  unvented  exploded  cavity  spheres  containing  a  fault  plane. 

The  configuration  for  simulating  a  vertical  fault  plane  is  siiown  in 
Figure  3.32(a).  The  6-inch-diameter  (15.2-cm)  fault  consists  of  three 
layers  of  low  friction  material  that  allow  for  slipping  along  the  plane. 
In  each  test,  the  fault  plane  passed  within  1-3/8  inch  (3.49  cm),  or 
approximately  1.8  exploded  cavity  diameters,  of  the  center  of  the  sphere. 

The  hydrofracture  record  for  test  192  is  shown  in  Figure  3.33.  The 
initially  low  cavity  pressure  indicates  that  charge  detonation  resulted 
in  venting  of  the  detonation  products.  The  subsequent  increase  in 
pressure  to  a  maximum  of  2660  psi  (18.34  MPa)  was  erratic  and  suggests 
that  fluid  was  escaping  from  the  cavity.  Examination  of  the  sphere 
following  the  test  revealed  that  the  crack  shown  in  Figure  3.34  extended 
from  the  cavity  to  the  fault  plane. 

The  hydrofracture  records  for  fault  plane  tests  210  and  221  are 
also  shown  in  Figure  3.33.  The  similarity  of  these  records  to  those  of 
tests  159  and  170  from  the  reproducibility  series  suggests  that  the  fault 
plane  was  not  a  factor  in  these  hydrofractures.  Cavity  gases  were  not 
vented  upon  charge  detonation.  Furthermore,  examination  of  each  sphere 
following  the  tests  revealed  that  fluid  escaped  from  the  side  of  the 
cavity  opposite  the  fault.  Fluid  then  flowed  around  the  annular  region 
of  residual  stress  and  eventually  reached  the  fault  later  in  the  hydro¬ 
fracture  process. 

Simulation  of  a  fault  plane  intersecting  a  cavity  was  provided  by 
the  configuration  shown  in  Figure  3.32(b).  The  3-inch-diameter  (7.62-em) 
fault  allows  for  slipping  along  the  plane.  The  resulting  hydrofracture 
record  for  test  217  is  shown  in  Figure  3.35.  For  comparison,  tests  159 
and  170  from  the  reproducibility  series  are  also  shown.  Initial  cavity 
gas  pressure  for  the  fault  test  is  in  the  range  of  pressures  for  tests 
without  the  fault.  Hence  cavity  gas  appeared  to  be  contained  as  hvdro- 
fracture  began.  As  hydrofracture  continued,  however,  cavity  pressure 
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FIGURE  3.32  CONFIGURATIONS  FOR  FAULT  PLANE  TESTS 
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170,  192,  210,  AND  221  —  NEARBY  VERTICAL  FAULT 


MP-8113-40 

FIGURE  3.34  HYDROFRACTURE  FROM  UNVENTED  EXPLODED  CAVITY 
TEST  192  NEARBY  VERTICAL  FAULT 
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FIGURE  3.35  HYDROFRACTURE  PRESSURES  FOR  UNVENTED  EXPLODED  CAVITY  TESTS  159, 
170,  AND  217  —  FAULT  THROUGH  CHARGE 


increased  slowly  and  erratically  to  a  plateau  at  2230  psi  (15.38  MPa), 
indicating  a  continual  loss  of  fluid.  The  sudden  loss  of  cavity  pressure 
at  the  end  of  the  plateau  corresponds  to  the  detection  of  a  surface 
crack.  The  absence  of  a  distinct  fracture  initiation  spike  in  the  pressure 
record  suggests  that  fluid  loss  occurred  along  the  fault.  Figure  3.36 
shows  that  a  fracture  extended  to  the  surface  of  the  sphere  along  the 
plane  of  the  fault. 
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Appendix  A 


FRACTURE  MECHANICS  ANALYSIS 


A  fracture  mechanics  analysis  is  presented  here  as  an  aid  to  inter¬ 
preting  the  hydrofracture  records  of  the  unexploded  cavity  tests.  The 
analytical  results  are  formulas  for  the  pressure  in  a  penny-shaped  crack 
that  will  cause  crack  propagation  in  an  infinite  elastic  medium.  The 
formulas  are  used  to  show  the  relationship  between  pressure  distribution 
and  crack  stability.  The  pressure  is  taken  as  uniform  over  a  constant 
central  circular  area  that  represents  the  unexploded  cavity,  and  different 
pressure  distributions  are  considered  in  the  surrounding  annular  region 
that  represents  a  crack  propagating  from  a  spherical  cavity.  Pressures 
required  to  start  expansion  of  the  circular  crack  from  the  constant  cen¬ 
tral  circular  region  are  modified  by  predictions  of  the  plate  edge  notch 
formulas  that  are  presented. 


The  stress  intensity  factor  for  a  penny-shaped  crack  of  radius  a, 
subjected  to  an  axisymmetric  pressure  distribution  p(r)  inside  the  crack 
is1! 


K 


I 


(A- I) 


A  concentric  region  of  radius  b  equal  to  that  of  the  spherical  cavity 
in  tiie  experiments  is  taken  to  be  loaded  by  a  uniform  pressure  distri¬ 
bution.  The  outer  annular  region  is  subjected  to  one  of  two  loadings: 
(1)  the  same  pressure  as  the  central  region  so  that  the  pressure  is 
uniform  throughout  the  crack  or  (2)  a  pressure  that  decreases  linearly 
from  the  central  region  pressure  to  zero  at  the  crack  tip.  Expressed 
mathematically,  these  two  cases  are 


A- 1 


. —  * 


p(r)  = 


P 


0  <  r  • 
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(A- 2) 


P(r)  = 


0  -  r  <  b 

b  <  r  <  a 


(A-3) 


Substituting  pressure  distributions  (A-2)  and  (A-3)  into  the 

integral  in  formula  (A-l)  and  integrating  provides  the  corresponding 

stress  intensity  factors.  As  the  pressure  P  is  increased  in  each  case, 

the  stress  intensity  factor  is  increased.  When  the  pressure  reaches  the 

critical  value  P  that  makes  the  stress  intensity  factor  equal  to  the 
e 

fracture  toughness,  K  ,  a  material  property,  the  crack  is  about  to  grow. 

If  an  incremental  increase  in  P  results  in  an  incremental  increase  in 

c 

the  crack  radius  a,  while  imposing  the  same  pressure  distribution,  the 
crack  growth  is  stable.  If  an  incrementally  larger  crack  radius  requires 
an  incrementally  smaller  pressure  to  maintain  =  K  ,  the  crack  growth 
is  unstable.  For  the  above  two  loading  cases  substitution  in  (A-l), 
integration,  and  use  of  gives  the  critical  pressures  P(  in  the 

f  orm 

-  L_ 

F^(a)  (case  1-uniform)  (A-4) 

—  (case  2- 1  inear )  (A-5) 

F2(-x) 

In  (A-4)  and  (A-5), 


F^(a)  =  a' 


(A-6) 


F2(a)  = 
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+  F  3 ( .  x )  (A- 7) 
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whore  a  =  a/b.  According  to  the  model,  b  is  a  constant  equal  to  the 
spherical  cavity  radius,  so  an  increase  of  a  from  its  initial  value  of 
unity  means  an  increase  of  the  crack  radius  a.  When  a  =  b,  that  is, 
when  x  =  1,  the  crack  has  a  radius  b  and  is  subjected  to  a  uniform 
pressure  over  its  entire  surface,  so  the  starting  critical  pressure  is 
given  by  the  classical  result 


x, 

k7  ("V  • 


(A-8) 


When  a  =  1,  the  functions  (A-b)  and(A-7)  each  have  the  value  unity,  which 
makes  formulas  (A-4)  and  (A-5)  initially  consistent  with  formula  (A-8). 

Figure  A.l  shows  the  variation  oi  critical  pressure  with  crack  radius. 
In  case  1,  a  uniform  pressure  always  acts  over  the  entire  crack  surface. 
Such  a  condition  may  be  realized  approximately  by  pressurizing  fluid  of 
very  low  viscosity  or  by  using  a  very  slow  pumping  rate.  The  crack  growth 
is  unstable  in  that  larger  cracks  require  less  pressure  just  to  keep  them 
stationary.  In  case  2,  the  pressure  applied  in  the  crack  decreases 
linearly  from  the  uniform  pressure  in  the  initial  penny-shape  to  zero 
at  the  crack  tip.  Such  a  condition  may  be  realized  approximately  by 
pressurizing  a  very  viscous  fluid.  The  crack  growth  is  initially  stable 
until  the  crack  radius  is  about  2-1/2  times  the  initial  radius.  Further 
growth  is  unstable. 

It  has  been  shown11  that  the  displacement  normal  to  the  face  of  a 
penny-shaped  crack  that  is  internally  pressurized  is 
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where  K  and  V  are  Young's  modulus  and  Poisson's  ratio  for  the  elastic- 
medium.  For  a  uniform  pressure  of  case  1,  formula  (A-9)  reduces  to 

(A-10) 


w(r) 
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v?)aP 
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which  gives  a  crack  volume  of 


V  =  16(1  -  v?)a3P/3E 


(A-l 1 ) 


The  critical  pressure-crack  radius  relationsliip  is  given  by  (A-4)  and 
(A-6),  so  elimination  of  the  radius  in  (A-ll)  leads  to  the  pressure- 
volume  relationship 
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This  relationship  is  shown  in  Figure  A. 2  for  =  400  psi  /in. 

(0.44  MPa  /m).  For  comparison,  the  water  hydrofracture  record  for 
test  165  conducted  at  zero  overburden  is  also  shown.  Theoretical  and 
experimental  data  both  indicate  a  monotonic  decay  of  cavity  pressure 
following  fracture  initiation.  Closer  agreement  between  the  curves 
would  be  realized  if  pore  fluid  migration  were  included  in  the  fracture 
mechanics  analysis. 

A  determination  of  fracture  initiation  pressure  requires  an  assumption 
on  the  initial  size  of  the  crack,  that  is,  the  initial  flaw  size  inherent 
in  the  material  or  caused  by  the  manufacturing  process.  If  the  flaw  is 
idealized  to  a  circular  notch  around  the  cavity  and  if  the  depth  of  the 
notch  is  small  relative  to  the  cavity  radius,  an  approximation  to  the 
stress  intensity  factor  may  be  taken  as  that  for  an  edge  notch.  The 
tensile  stress  acting  on  the  notch  is  the  circumferential  stress  at  the 
cavity  caused  by  internal  and  external  pressures,  p  and  p^.  As  shown 
in  Figure  A. 3(a),  this  stress  is  p  ^ / 2  -  3p^/2,  and  the  stress  intensity 
factor  is,? 
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FIGURE  A  2  PRESSURE-VOLUME  RELATIONSHIP  FOR  A  UNIFORMLY  PRESSURIZED  CRACK 


K]  =  1.12  Pa  -  |  pbJ  [  ii  ( a  -  b)]  '2  (A-13) 

For  a  constant  fluid  pressure  p  in  an  edge  notch,  as  shown  in 

a 

Figure  A.  3(b)  the  stress  intensity  factor  is1*’ 


Kr  =  1.12  pa  [ it ( a  -  b)]  '2 


(A— 14) 


Superposition  of  the  loading  states,  as  shown  in  Figure  A. 3(c)  gives  the 
stress  intensity  factor  that  is  the  sum  of  (A-13)  and  (A-14),  namely. 
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The  critical  pressure  difference  is  therefore,  from  (A— 15), 
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The  above  analysis  is  now  used  to  explain  the  dependence  of  fracture 
initiation  pressure  on  hydrofracture  fluid  viscosity  and  overburden  as 
shown  in  Figure  3.5  of  the  main  text  of  this  report. 

All  tests  were  performed  with  a  constant  flow  rate  of  4.26  cm 5 /min 
into  3/4-inch-diameter  (1.90-cm)  smooth  cavities.  The  fracture  initiation 
pressures  for  glycerol  are  well  above  those  for  water.  We  postulate  that 
at  the  cavity  wall  an  edge  notch  exists  that  is  subjected  to  a  uniform 
internal  pressure,  which  is  the  pressure  in  the  hydrofracture  fluid.  If 
the  overburden  and  cavity  pressures  are  p  and  p  and  the  notch  length 
is  c,  the  stress  intensity  factor  Kj  is 

Kj  =  1.68  (p.(  -  pb)(nc)  2  (A- 1  7 ) 


A-8 


From  the  results  for  water  in  Figure  3.5,  p  -  p  =  800  psi  (5.52  MPa) 

ti  D 

and  from  fracture  tests  on  2C4  the  fracture  toughness  is  =  400  psi 
/  in .  (0.44  MPa/in").  If  we  set  Kj  =  K ^  and  use  the  experimentally 

determined  pressure  difference,  equation  (A-17)  gives  a  constant  crack 
length  of  c  =  28  mils  (0.71  mm). 

There  is  no  reason  to  assume  that  the  inherent  initial  crack  length 
is  different  in  the  2C4  grout  spheres  that  were  hydrof ractured  with 
glycerol.  The  hydrostatic  fracture  mechanism  just  outlined  leads  to 
the  same  fracture  initiation  pressure  in  glycerol  as  that  in  water.  It 
is  postulated,  therefore,  that  glycerol  initiates  fracture  at  the  same 
uniform  pressure  as  water,  but  immediately  after  the  crack  tip  moves 
and  allows  the  glycerol  to  commence  flow  in  the  crack,  the  viscosity 
causes  a  substantial  redistribution  of  pressure  in  the  crack.  The 
pressure  falls  to  a  low  value  near  the  crack  tip  and,  because  of  the 
constant  pumping  rate,  it  rises  in  the  cavity.  This  pressure  redistri¬ 
bution  creates  a  lower  stress  intensity  factor.  As  an  example,  we 
approximate  the  pressure  distribution  by  a  linear  decrease  from  p  in 
the  cavity  to  p^  near  the  crack  tip.  The  stress  intensity  factor  is 
then 

Kt  =  [1.68  (p  -  p  )  -  0.68  (p  -  P  )  ]  (ttc )  (A- 18) 

1  ci  D  <J  C 

Assuming  that  the  crack  has  moved  only  a  small  amount  to  start  flow 
along  the  crack,  we  see  that  when  the  experimental  values  of  p^  =  2700  p 
(18.62  MPa)  and  p^  =  1000  psi  (6.90  MPa)  again  initiate  crack  growth  we 
have  p^  =  483  psi  (3.33  MPa)  near  the  crack  tip.  If  the  crack  wore 
49  mils  (1.24  mm)  long  at  this  stage  instead  of  28  mils  (0.71  mm)  long, 
we  would  have  p^  =  0  near  the  crack  tip  when  crack  growth  was  restarted 
by  =  2700  psi  and  =  1000  psi.  For  radial  laminar  flow  in  a  crack 
of  constant  thickness  h  and  length  c  encircling  a  cavity  of  radius  a, 
the  pressure  drop  from  the  cavity  to  the  crack  tip  is 

p  =  12p(Q/2fta)c/h  '  (A- 19) 


A-9 


Hence  a  pressure  drop  of  2700  psi  (18.62  MPa)  in  a  crack  of 
49  mils  (1.24  mm)  extending  from  a  cavity  of  3/8  inch  (0.95 
implies  a  crack  thickness  of  h  =  0.3  mils  (0.0076  mm).  All 
numbers  in  the  example  seem  practical. 


1 ength 
cm)  radius 
the  above 


If  the  crack  in  the  cavity  is  approximated  by  a  penny-shaped  crack, 
a  pressure  distribution  in  the  crack  that  decreases  from  the  cavity 
pressure  to  zero  at  the  tip  produces  stable  crack  growth  followed  by 
unstable  crack  growth  as  shown  in  Figure  A.l;  thus,  some  growth  is 
possible  during  the  initial  pressure  rise  of  the  hydrofracture  record. 
Such  stable  growth  probably  accounts  for  the  maximum  slope  of  the  initial 
pressure  rise  portions  of  the  hydrofracture  record  frequently  being  at 
pressures  below  the  peak  pressure. 

Finally,  differentiation  of  (A-17)  yields  the  theoretical  slope  of 
the  fracture  initiation  versus  overburden  curve  for  both  water  and 
glycerol  hydrofracture,  that  is. 
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which  is  in  agreement  with  the  experimental  data  for  RMG  2C4  shown 
earlier  in  Figure  3.5. 
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Appendix  B 


POROUS  FLOW  ANALYSIS 

A  porous  flow  analysis  based  on  the  work  of  Rice  and  Cleary1 j  is 
presented  here  as  an  aid  to  understanding  the  decay  of  an  explosively 
generated  residual  stress  field  surrounding  a  spherical  cavity.  The 
analytic  results  are  formulas  for  the  stress-induced  flow  of  interstitial 
fluid  in  a  porous  solid.  The  formulas  are  used  to  show  the  relation¬ 
ships  between  pore  pressure  and  the  effective  stresses  acting  on  the 
solid  matrix.  The  initial  pore  pressure  distribution  is  approximated 
by  the  mean  normal  stress  associated  with  a  calculated  residual  stress 
field.  Equations  are  specialized  for  the  case  of  spherically  symmetric 
mot  ion . 

The  total  stresses  ck.  acting  on  a  material  element  and  the  pore 
pressure  p  are  defined  as  the  basic  state  variables  and  are  assumed 
related  in  some  appropriate  fashion  to  the  strains  of  the  solid  and 

the  mass  m  of  pore  fluid  per  unit  reference  volume.  The  linearized 
constitutive  relations  for  a  fluid-saturated  porous  elastic  solid  are 
presented  for  quasi-static  elastic  deformation  in  which  local  equilibrium 
of  the  pore  fluid  can  be  assumed. 

For  the  strains  of  the  solid 
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where  l>  and  v  are  the  shear  modulus  and  Poisson  ratio  when  the  material 
is  deformed  under  drained  conditions,  v  is  the  undrained  Poisson  ratio, 
and  B  is  defined  by  the  relation  between  initial  induced  pore  pressure 
and  the  total  hydrostatic  stress  acting  on  an  element.  That  is. 
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The  elastic  constants  v  and  B  may  be'  determined  directly  from  the  result 
of  a  single  undrained  test  in  which  the  Poisson  effect  and  induced  pore 
pressure  are  measured. 

For  the  mass  of  pore  fluid  per  unit  reference  volume 
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where  m  and  c  are  the  fluid  mass  content  of  an  element  and  the  fluid 
o  o 

density  for  the  unstressed  state. 


is 


The  constitutive  law  governing  pore  fluid  diffusion,  D'Arey's  law, 
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where  q  is  the  fluid  mass  flow  rati1  in  the  radial  direction  per  unit 
area  and  k  is  the  permeability. 

The  field  equations  are  also  developed  for  quasi-static-  response. 
Kqui librium  of  stresses  requires  that 
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The  kinematic-  relations  between  strain  and  displacement  u  in  the  radial 
direct  ion. 
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lead  to  the  compatibility  equations  as  follows: 
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Use  of  equilibrium  (B-6)  and  the  constitutive  relations  (B-l  and  B— 2 )  in 
the  compatibility  condition  (B-8)  yields 
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where 
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Hence  from  (B-9) 
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Tile  final  governing  equation  is  for  mass  conservation  of  the  infiltrating 
pore  fluid 


3r 


a 

+  - 
r 


q 


3m 

3t 


(B-l 2) 


B-l 


Is.  of  C H---+ )  and  (B-5)  in  (B-12)  yields 
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gives 

/  ;V p  2  dp  N 

1  y  dr-  r  dp 

\  =  Ip.  +  — "''ii 

/  dt  (1  -  V  )(1 

U 

,')  dCj(t) 

+  v )  d  t 

the  coefficient  of  consolidation  and  is  given  by 
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I’sl*  of  (B-ll)  to  eliminate  in  (B-6)  ,  and  integration  of  the  result 

subject  to  the  condition  c  (r  =  a)  =  o  (a)  vields 
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liquation  (B-14)  governs  pore  pressure  distribution.  liquations  ( B  —  1  1 ) 
and  (B-lb)  govern  the  change  in  stress  associated  with  a  change  in  pore 
pressure.  For  the  special  case  of  an  infinite  medium,  o  (»)  =  '  (■)  = 
p(')  =  0.  Hence  (B-ll)  requires  that  C^(t)  =  0. 

We  apply  these  results  assuming  that  the  2C4  properties  are  not 
changed  by  the  explosion.  Consider  the  infinite  region  surrounding  an 
exploded  cavity.  Assume  that  the  step  function  shown  in  Figure  B.l(a) 
approximates  the  initial  pore  pressure  resulting  from  charge  detonation. 

The  height  of  the  step  is  determined  from  the  mean  normal  stress  associated 
with  a  calculated  residual  stress  field.'1  The  width  of  the  step  is  taken 
as  the  extent  of  the  undyed  region  associated  with  hvdrof racture  of  an 
exploded  cavity  sphere.  Hence  the  initial  pore  pressure  distribution  is 
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The  radial  and  tangential  stresses  associated  with  (B-17)  are  determined 
from  (B-ll)  and  (B-lb)  with  C^(t)  =  0: 
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The  pore  pressure  distribution  for  t  '  0  is  given  by  the  solution  (see 
Reference  14)  of  (B-14)  with  C^(t)  =  0: 
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2r/  net 


—  ( r— ( >)  ‘  /4ct 


- ( r+,  -2a) ' / 4c  t 


( B-22) 


For  times  long  enough  sucli  that 


2 

4  c.t  >>  (b  -  a) 


B-6 


(B-2-1) 


p 


the  following  approximations  are  valid: 


e-(r-p)-'/4ct 


(r  -  ,  )  ‘  /Act 


and 


-(r+4  -2a)  ’  /  Ai  t 


(r  +  .  -  2a)  '/Aet 


(B-2A) 


(B-25) 


Substituting  (B-2A)  and  (B-25)  into  (B-22)  and  integrating  yields  the 
following  approximation  for  late-time  pore  pressure  distribution: 


P(r,tco)  » 


IPo_  (b  -  nj 
3/n  Act 


(b  a) 
►  Act 


2  b  +_  a 
b  -  a 


(B-26) 


Numerical  values  of  relevant  material  properties  are  shown  in 
Table  B.l.  Drained  and  undrained  shear  modulus  and  Poisson  ratio  are  based 
on  acoustic  measurements  performed  on  a  saturated  and  partially  dried 
specimen  of  RMG  2CA.  Permeability  was  measured  in  a  water  permeameter.  Ha- 
value  for  B  was  estimated  from  a  table  of  geological  properties  presented 
in  Rice  and  Clary.1  '  The  coefficient  of  consolidation  and  n  were  determined 
from  ( B— 15)  and  (B-10),  respectively. 


Table  B.  1 

MATERIAL  PROPERTIES  FOR  POROUS  FLOW  ANALYSIS 


Property _ RMG  2CA 


Drained  shear  modulus  G 

(psi ) 

1.08  x 

101 

Drained  Poisson  ratio  v 

0.275 

Undrained  shear  modulus 

G  (psi) 

U 

1.03  x 

10 

Undrained  Poisson  ratio 

\1 

u 

0.282 

Permeability  K(em,t /dyne- 

sec  ) 

1  x  10" 

1 

B 

0.6 

Coefficient  of  consolidation  c(cm’/sec) 

A.2A 

n 


0.019 


The  region  of  significant  initial  pore  pressure  is  specified  by 
taking  a  =  1  cm  and  b  =  4  cm  in  (B-17).  From  (B-23),  the  late- time 
approximation  is  valid  for 


4ct 

(b  -  a)2 


10 


(15-27  ) 


or 


t  =5.3  sec 

CO 


(B-28) 


From  (B-26),  the  corresponding  late-time  pore  pressure  distribution  is 
given  by 

P(r,tJ  -  0.036  po(l  -  *)  (B-29) 

and  is  shown  schematically  in  Figure  B.l(b).  If  the  total  residual 
stress  field  generated  by  charge  detonation  is  approximated  by 

=  p^,  then  the  maximum  percentage  decay  due  to  pore  fluid 
migration  is  found  from  (B-18)  and  (B-20)  as  follows: 


and 


l1r  (b,0) 


On(b,0) 


4,1 

3 


2n 

3 


(■ 

(' 


3 

a 


) 

) 


=  2.5% 


=  1.2% 


(13-30) 


(B-31 ) 
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Appendix  C 


MATERIAL  PROPERTIES 

Table  C.l  provides  the  mixtures  and  heat  curing  cycles  for  rock¬ 
matching  grout  RMG  2C4,  low  density  rock-matching  grout  LD  2C4,  and 
granite  simulant  GS3.  Table  C.2  summarizes  the  commonly  used  physical 
and  mechanical  properties. 

The  scatter  in  crush  strength,  and  tensile  strength  associated  with 
the  testing  of  a  large  number  of  2-inch-diameter  (5.08-cm)  cylinders  is 
shown  in  Figures  C.l  through  C.4.  These  figures  show  that,  in  general, 
reasonable  tolerance  is  maintained  on  material  strength.  Additional 
strength  properties  of  2C4  and  GS3  are  shown  in  Figures  C.3  and  C.6. 
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Table  C.l 


MIXTURES3  FOR  ROCK-MATCHING  GROUT  RMG  2C4, 
LOW  DENSITY  ROCK-MATCHING  GROUT  LD  2C4, 
AND  GRANITE  SIMU1.ANT  GS3 


Component 

RMC  2C4 

LD  2C4 

Type  I— I I  Portland  cement 

32.691% 

32.063% 

Sand  (20-40  Monterey) 

21.896% 

21.476% 

Barite  (barium  sulfate) 

20.848 % 

20.447% 

Bentonite  (gel) 

2.837% 

2.783% 

CFR  2  (concrete  friction 
reducing  compound) 

0.078% 

0.077% 

Water 

21.650% 

21.234% 

Microballoons  (Q-CEL  300) 

— 

1.920% 

Component 

GS3 

Type  I — I I  Portland  cement 

38.573% 

Sand  (granite) 

49.295% 

Me  Intent 

1  .035% 

CFR  2 (concrete  friction 
reducing  compound) 

0.252% 

Water 

10.845% 

‘'28-day  aging  procedure:  sealed  in  plastic  and  submerged 
in  water  with  the  following  temperature  sequence:  raise 
to  54°C  over  48-hour  period,  hold  at  54°C  for  48  hours, 
lower  to  25°C  over  36-hour  period. 


C-2 


Table  C.2 


PROPERTIES  OF  ROCK-MATCHING  GROUT  RMG  2C4 , 
LOW  DENSITY  ROCK-HATCHING  GROUT  LD  2C4, 
AND  GRANITE  SIMULANT  GS3 


(a)  PHYSICAL  PROPERTIES 


Physical  Property 

RMG  2C4 

LD  2C4 

GS3 

3 

Density  (g/cm  ) 

Aged 

2.15 

1.90 

2.42 

Dry 

1.75 

1.57 

2.27 

Grain 

2.87 

2.68 

2.80 

Water  by  wet  weight  (%) 

18.6 

17.4 

6.3 

Porosity  (%) 

39 

43 

18.8 

Saturation  (%) 

100 

86 

81.6 

Air  voids  (%) 

0 

13.4 

3.5 

Longitudinal  velocity  (km/sec) 

3.29 

3.13 

4.82 

Shear  velocity  (km/sec) 

1.82 

1.78 

2.75 

Modulus  in  compression  (psi) 

2.64  x  106 

2.20  x  106 

6.68  x  106 

Shear  modulus  (psi) 

1.03  x  106 

8.74  x  105 

2.65  x  106 

Bulk  modulus  (psi) 

2.00  x  106 

1.53  x  ]06 

4.64  x  106 

Poisson's  ratio 

0.28 

0.26 

0.26 

Permeability  (pd) 

3.0 

- 

- 

(b)  MECHANICAL  PROPERTIES 


Material 

Average 

Strain  Rate 
(sec-l) 

Compressive 

Strength 

(psi) 

Tens  i  1  e 
Strength 
(psi) 

J- 

RMG  2C4 

Static 

3970 

530 

RMG  2C4 

0.15 

5330 

900 

LD  2C4 

Static 

3200 

460 

LD  2C4 

0.15 

5000 

780 

OS  3 

Static 

8900 

1110 

GS3 

0.15 

1  5900 

1620 
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FIGURE  C.1  UNCONFINED  CRUSH  STRENGTH  OF  ROCK-MATCHING  GROUT 
RMG  2C4 


C-4 


-•->*  •  * vr_  - 


n 


SPLITTING  TENSILE  STRENGTH  —  MPa 


0  1  2  3  4  5 
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MA-81 13-50 

FIGURE  C.2  SPLITTING  TENSILE  STRENGTH  OF  ROCK-MATCHING  GROUT 
RMG  2C4 


C-5 


NUMBER  OF  TESTS 


UNCONFINED  CRUSH  STRENGTH  —  MPa 
0  10  20  30  40  50  60  70 


MA-8113-51 


FIGURE  C.3  UNCONFINED  CRUSH  STRENGTH  OF  GRANITE  SIMULANT  GS3 
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FIGURE  C.4  SPLITTING  TENSILE  STRENGTH  OF  GRANITE  SIMULANT  GS3 
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STRESS  DIFFERENCE,  a o3  —  kbar 


0  0.1  0.2  0.3  0.4  0.5  0.6 


CONFINING  PRESSURE,  o3  —  kbar 
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FIGURE  C.5  STRENGTH  PROPERTIES  OF  ROCK-MATCHING  GROUT  RMG  2C4 
AND  GRANITE  SIMULANT  GS3:  STRESS  DIFFERENCE  VERSUS 
CONFINING  PRESSURE 
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Appendix  D 


SEQUENCE  OF  EVENTS  IN  UNVENTED  EXPLODED  CAVITY  TESTS 

Figure  D.l  shows  a  schematic  representation  of  cavity  pressure  and 
fluid  line  pressure  generated  during  a  typical  unvented  exploded  cavity 
tost.  Charge  detonation  occurs  at  time  zero  and  results  in  an  initial 
cavity  gas  pressure  indicated  by  point  A.  Curve  AB  represents  the  decay 
of  residual  gas  pressure  in  a  closed  cavity  as  the  detonation  products 
cool  and  diffuse  into  the  surrounding  medium.  Point  C  represents  the 
prepressure  applied  to  the  hydrofracture  fluid.  Pumping  begins  at  time 
zero,  and  the  straight  line  segment  CD  represents  pressurization  of  the 
fluid  line  as  pumping  continues.  The  ball  valve  that  seals  the  cavity 
starts  to  open  at  time  zero.  However,  cavity  gas  pressure  is  high  enough 
for  four  seconds  following  charge  detonation  to  prevent  the  valve  from 
opening.  Curves  BE  and  DE  represent  equilibration  of  cavity  gas  pressure 
and  fluid  line  pressure  when  the  cavity  finally  opens.  Hydrofracture 
records  for  unvented  exploded  cavity  tests  ,->tart  at  point  E,  and  curve  EF 
represents  the  initial  portion  of  a  typical  record. 

The  300  psi  (2.068  MPa)  at  point  C  is  the  nominal  prepressure  applied 

to  the  hydrofracture  fluid  before  charge  detonation.  The  1176  psi 

(8.108  MPa)  at  point  D  is  the  pressure  generated  in  the  fluid  line  after 

four  seconds  of  pumping  at  a  rate  of  4.26  cmVmin.  This  value  is  obtained 

3 

from  a  calibration  test  on  the  closed  system  of  90  cm  of  water.  The 

3 

pressure  in  the  closed  system  increases  3089  psi  (21.30  MPa)  for  every  cm 

of  fluid  compression.  The  1500  psi  (10.34  MPa)  at  point  E  is  the  measured 

value  of  equilibrium  pressure  between  cavity  gases  and  hvdrofracture  fluid 

when  the  cavity  first  opens.  The  1550  psi  (10.69  MPa)  at  point  B  is  a 

theoretical  estimate  of  cavity  gas  pressure  just  before  the  cavity  opens. 

The  calculation  is  based  on  the  assumption  that,  when  the  cavity  opens, 

3 

the  3.6  cm  of  detonation  products  in  the  cavity  expand  as  a  perfect  gas 
with  a  ratio  of  specific  heats  of  1.4.  ‘ 
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FIGURE  0.1  SEQUENCE  OF  EVENTS  DURING  UNVENTED  EXPLODED  CAVITY  TESTS 


r 


i 


Tlu'  liydrof racture  records  for  unvented  exploded  cavity  tests  indicate 

that  the  pressure  at  fracture  initiation  depends  on  the  minimum  residual 

cavity  gas  pressure.  As  an  aid  to  interpreting  this  dependence,  Figure  D.2 

shows  the  stress  history  at  the  cavity  wall  together  with  a  Mohr-Coulomb 

yield  surface.  Charge  detonation  results  in  a  cavity  pressure  increase 

indicated  by  the  path  0-1-2.  Point  2  corresponds  to  point  A  in  Figure  D.l. 

Cavity  gases  then  cool  along  path  2-3-41  Point  4'  corresponds  to  the 

minimum  cavity  pressure.  Hydrofracture  results  in  a  cavity  pressure 

increase  along  path  4  *  —5 * .  Cavity  fracture  occurs  at  point  5'  where  tlu* 

circumferential  stress  a„  equals  the  static  tensile  strength  0„  of  the 

0  1 

material .  Fracture  does  not  occur  along  tiie  path  0-1-2  because  the  dynamic 
tensile  strength  is  large  enough  to  withstand  the  dynamic  loads.  Figure  D. 
shows  that,  if  the  cavity  pressure  had  been  unloaded  to  point  4"  instead 
of  point  4',  the  resulting  fracture  initiation  pressure  would  have  been 
given  by  point  5"  instead  of  point  5'.  That  is,  a  lower  residual  cavity 
gas  pressure  results  in  a  lower  fracture  initiation  pressure. 
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FIGURE  D.2  STRESSES  AT  CAVITY  WALL  DURING  UNVENTED 
EXPLODED  CAVITY  TESTS 
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